The 


Canadian Journal 
of 


Chemical Engineering 


formerly 


CANADIAN JOURNAL OF TECHNOLOGY 


CONTENTS 


The Upwards Vertical Flow of Air-Water Mixtures G. W. Govier 
Il. Effect of Tubing Diameter on Flow Pattern, W. Leigh Short 


Holdup and Pressure Drop 


A Heat Exchanger Design for Heavy Water Reactor A. D. Duff, Jr. 
Service E. E. Wilson 


First Order Rate Processes and Axial Dispersion in J. J. Carberry 
Packed Bed Reactors 


Correction Factor for Axial Mixing in Packed Beds N. Epstein 


Some Procedures for the Evaluation of Reactor Fuels R. F. S. Robertson 
and Sheathing Materials at Chalk River 


Conditioning of D,O in Heavy Water Power Reactors G. M. Allison 


End Effect Corrections in Heat and Mass Transfer A. 1. Johnson 
Studies A. Hanzielec 
D. Ward 
A. Golding 


Viscous Flow in Multiparticle Systems John Happel 








A Pigmy Paper-maker 


turns out samples at the C-I-L Central Research Laboratory 










at Beloeil, Que. The products of this miniature paper mill 


. : : ; : | 
are used to study bleaching of various types of wood pulp with ) 
C-I-L chemicals. The findings of the studies are made | 


available to the pulp and paper industry through the technical 


service representatives of C-I-L. | 


C-I-L research helps many other industries to develop the 


improved products Canada needs for progress. 


CANADIAN INDUSTRIES LIMITED | 


Serving Canadians Through Chemistry | 





Agricultural Chemicals * Ammunition * Coated Fabrics ~- Industrial Chemicals * Commercial Explosives * Paints « Plastics * Textile Fibres 





Fibres 


The Canadian Journal of Chemical Engineering 


Jormerly 


Canadian Journal of Technology 





VOLUME 36 


OCTOBER, 1958 


NUMBER 5 





EDITOR: 


W. M. CAMPBELL 


Chemistry and Metallurgy Division, Atomic Energy of Canada Limited, 
Chalk River, Ont. 


Managing Editor: 
T. H. G. MICHAEL 


Editorial Assistant: 


R. G. WATSON 


Assistant Editor: 
R. N. CALLAGHAN 


Circulation Manager: 


M. M. HOLDEN 


EDITORIAL BOARD: 


Chairman 
J. R. DONALD, Montreal, Que. 


A. CHOLETTE, Quebec, Que. 
W. H. GAUVIN, Montreal, Que. 
GLEN GAY, Ottawa, Ont. 

G. W. GOVIER, Edmonton, Alta. 


E. B. LUSBY, Toronto, Ont. 
LEO MARION, Ottawa, Ont. 
R. R. McLAUGHLIN, Toronto, Ont. 


G. L. OSBERG, Ottawa, Ont. 
E. R. ROWZEE, Sarnia, Ont. 





| A. I. JOHNSON, Toronto, Ont. 
H. R. L. STREIGHT, Montreal, Que. 


EX-OFFICIO: 


H. BORDEN MARSHALL, Toronto, Ont. 
W. G. DICKS, Beloeil, Que. 


C. E. CARSON, Toronto, Ont. 
H. S. SUTHERLAND, Montreal, Que. 





Authorized as second class mail, Post Office Department, Ottawa. Printed in Canada 


| 


Manuscripts for publication should be submitted to the 
Editor: Dr. W. M. Campbell, Chemistry and Metallurgy 
Division, Atomic Energy of Canada Limited, Chalk River, 
Ontario. (Instruction to authors are on inside back cover) 


Proofs, correspondence concerning proofs, and orders 
for reprints should be sent to: The Chemical Institute of 
Canada, 18 Rideau Street, Ottawa 2, Ont. 


Subscriptions, renewals, requests for single or back 
numbers and all remittances should be sent to: The Chemical 
Institute of Canada, 18 Rideau Street, Ottawa 2, Ont. 


C.J.Ch.E. is published by The Chemical Institute of 
Canada every two months. Subscription rate is $3.00 per year 
and .75¢ per single copy; U.S. and U.K.—$4.00; and Foreign— 
$4.50 per year. 


The Canadian Journal of Chemical Engineering, October, 1958 


Unless it is specifically stated to the contrary, the Institute 
assumes no responsibility for the statements and opinions 


expressed in The Canadian Journal of Chemical Engineering. 
Views expressed in the editorials do not necessarily represent 
the official position of the Institute. 


Editorial Offices: Dr. W. M. Campbell, Chemistry and 
Metallurgy Division, Atomic Energy of Canada Limited, 
Chalk River, Ont. 


Production and circulation offices: 18 Rideau Street, 


Ottawa 2, Ont. 


Change of Address: Advise the Circulation Department, 
The Chemical Institute of Canada, 18 Rideau Street, Ortawa 2, 
Ont. in advance of any change of address, providing old as 
well as new address. Enclose address label if possible. 


193 

















The Upward 






Vertical Flow of Air-Water Mixtures’ 


II. Effect of Tubing Diameter on Flow 
Pattern, Holdup and Pressure Drop 


G. W. GOVIER? and W. LEIGH SHORT: 


Results of visual observations and measurements 
of holdup and pressure drop are reported over a 
range of air-water ratios for the upwards vertical 
flow of air-water mixtures in four tubes with diam- 
eters ranging from 0.630 in. to 2.50 in. Tests were 
conducted under conditions of constant air density 
corresponding with an average pressure of 36 psia 
and at 10 constant superficial water velocities ranging 
from 0.0695 to 7.35 ft./sec. 

The data are analyzed by the method previously 
suggested by Govier, Radford and Dunn involving the 
separation of the unit pressure drop curves into four 
regimes to aid in flow pattern description and the 
separation of the unit pressure drop itself into two 
components to facilitate correlation. 

Tubing diameter is found to have little or no 
effect on the border between Regimes I and II or 
on the transition from the bubble to the slug flow 
pattern. The borders between Regimes II and III and 
between Regimes III and IV are markedly affected 
by tubing diameter as by consequence are the transi- 
tions from slug to froth and from froth to ripple flow. 

Tubing diameter is found to have an important 
effect on the superficial friction factor and the hold- 
up. The pressure drop data for each of the diameters 
have been successfully brought together in a single 
friction factor-Reynolds number correlation. The 
holdup data are separately correlated for pressure 
drop Regimes I and II and for Regime III and Regime 
IV. 


2 first paper “) of this series presented a review of 
the general characteristics of the upwards vertical 
flow of gas-liquid mixtures and indicated the complex 
flow patterns and the existence of holdup both of which 
distinguish the flow of two phase mixtures from that of 
single phase fluids. A comparison of the flow pattern 
nomenclature employed by other workers led to the 
proposal of the terms “bubble”, “slug”, “froth”, “ripple”, 
“film” and “mist” for the identifiable patterns encountered 
at constant liquid flow rate with increasing ratio of the 
gas to liquid phase. These flow patterns were identified 
with some precision and correlated with the aid of the 
loci of inflection points in the relationships between the 
unit pressure drop and the volume ratio of the phases. 
Holdup (which is related to “slip velocity”) was ex- 
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pressed quantitatively in terms of holdup ratio, Ha, 
defined as the ratio between (a) the gas-liquid volume 
ratio in the supply mixture and (b) the gas-liquid volume 
ratio in the flow section. This holdup ratio was correlated 
in terms of the superficial velocity of the liquid phase 
and the volume ratio of the phases. 

In a thermodynamic analysis of the problem the overall 
unit pressure drop was separated into its two basic com- 
ponents—that due to the hydrostatic head and that due to 
the irreversibilities attending flow. The equation expres- 
sing this relationship, in terms of feet of liquid phase per 
foot of height, is: 


~y, AP _ 1+ Re ines i (AF 1 
ae "ie 14h 
where 
Vv, = specific volume of liquid, ft®/Ib 
—AP = pressure drop, lbs/ft? 
AX = incremental height, ft 
Ry = i ate mass ratio, ? 
Ry = gas-liquid volume ratio, pve 
L 
(Ae _ irreversibility effect, feet of liquid per 
AX /t foot of height 


The quantity (=) may conveniently be expressed in 
« L 


terms of f’;, a superficial friction factor based upon the 
liquid phase, by the relationship: 


AF 2 f’L V2t 
<i re oe ee (2) 
AX], g- D 
where 
V,, = the superficial liquid phase velocity based upon 
the total tube cross section, ft./sec. 
g. = dimensional conversion factor, lb..ft./Ib., sec.? 


D = tube diameter, ft. 

This latter equation serves to define the superficial friction 
factor, f’;, which may be expected to depend upon a 
superficial Reynolds number, and some function of the 
system geometry, the fluid properties and the mass or 
volume ratio of the two phases. The work previously 
reported has indicated that experimental pressure data 
may in fact be so correlated. 

The correlations of the flow pattern, holdup and pres- 
sure drop data referred to, however, are limited in appli- 
cations to the conditions of the experimental data. The 
data presented in Part I are for the air-water system 
flowing in a smooth bore 1.025 in. I.D. tube at water rates 
ranging from 0.00040 to 0.0421 cu. ft./sec. and at air-water 
volume ratios ranging from 0 to 348. The properties of 
both phases were held constant through the maintenance 
of constant average flowing pressure and temperature of 


195 











36.0 psia and 85 F. The work reported here was under- 
taken to shed light on the effect of tubing diameter. Other 
work now under way and to be reported later has been 
designed to clarify the role of liquid and gas phase pro- 
perties. 


Experimental equipment and procedure 

The equipment employed for the tests reported here 
and illustrated in Figure 1 is essentially the same as that 
described in Part I. It was modified by the incorporation 
of three additional flow sections to give the following 
each of 22.88 ft. effective length: 

(a) a 0.630 in. I.D. smooth bore transparent cellulose- 

acetate-butyrate tube, 

(b) the original 1.025 in. 1.D. smooth bore transparent 

cellulose-acetate-butyrate tube, 

(c) a 1.50 in. LD. smooth bore seamless copper tube, 

and 

(d) a 2.50 in. LL.D. smooth bore opaque carbon black 

loaded polythene tube. 
Suitable carefully reamed pressure taps were made in the 
new tubes at elevations corresponding to that in the 
original tube. 

The data taking procedure was almost exactly as pre- 
viously reported, but flow patterns were identified pri- 
marily from the inflection points. Tests were conducted 
at superficial water velocities at and near 0.0695, 0.135, 
0.295, 0.625, 0.869, 1.50, 2.05, 2.51 and 3.47 ft./sec. corre- 
sponding to those previously made in the 1.025 in. tube. 
The average flowing pressure was held constant at 36.0 
psia and the temperature was held nearly constant in the 
range 70°-85°F. 


To Pressure 
Controller 









Plug Valves a | 


Downstream Pressure 


Stations Air -Water 


Seporator 


Mid-Point Pressure Stations 
(To Pressure Controller) 
Water Rotameters 


i 


= 
\" Water Supply 
210 PS.LG. 






I" Air Supply 
75 PSIA.G. 





Filter 
Stations 


Plug Valves 


Air Orifice 
Metering Station 


Drain 


Figure 1—Schematic flow diagram of equipment. 


Experimental results 

The basic experimental data are tabulated in Table | 
(Partial) for the typical constant water rate of some 0.87 
ft./sec. The complete tabular data for all water rates is 
given in Table 1 deposited as Document No. 5723* with 
the American Documentation Institute Auxiliary Publica- 
tion Project, Photo Duplication Service, Library of Con- 
gress, Washington 25, D.C. For convenience the data from 
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Part I for the 1.025 in. tube at the 0.873 ft./sec. velocity 
is repeated in Table I (Partial). 

Referring to Table I (Partial), column 2 gives the 
tube diameter, column 3 gives the water rate in cubic feet 
per second, column 4 the water-saturated air rate in cu. 
ft./sec. at the mid-point temperature given by column § 
and the constant mid-point pressure of 36.0 psia. These 
rates were calculated from the actual measured air rates 
to correct for the saturation of the air with water which 
took place in the calming section. Column 6 presents the 
ratio of the air to water rates on a volume basis. Column 
7 gives the volume ratio of air to water trapped in the test 
section at the completion of each test. Column 8 gives 
the specific volume of the saturated air at mid-point test 
section conditions. Column 9 presents the holdup ratio, 
being the ratio of the discharge air-water volume ratio to 
the test section air-water volume ratio i.e. of column 6 to 
column 7. Column 10 gives the total pressure drop be- 
tween the upstream and downstream pressure measuring 
stations, expressed in terms of ft. of water per ft. of 
height. Column 11 identifies the pressure drop regime and 
Column 12 the flow pattern. 

The basic data for the 0.630, the 1.50 and the 2.50 in. 
I.D. tubes at each of the water rates tested are summarized 
in Figures 2 through 18. These figures show by separate 
curves (a) the relationship between the unit pressure drop 
and the air-water volume ratio at test conditions and (b) 
the relationship between holdup ratio and air water 
volume ratio at test conditions. The four pressure drop 
regimes defined in Part I for the 1.025 in. tube are clearly 
distinguishable in Figures 2, 3, 4 and 12, perceptible in 
Figures 5, 6 and 11, and partially apparent in Figures 7 to 
10, and 13 to 18. As previously found the data for the 
higher water rates indicates a narrowing of Regime III 
and a merging of Regimes II and IV. Thus for the 0.630 
in. tube all four regimes are apparent at a superficial water 
velocity of 0.624 ft./sec., Regime III is barely perceptible 
at the velocity of 2.04 ft./sec. and it has disappeared 
completely and Regimes II and IV have merged at the 
velocity of 3.47 ft./sec. Similar observations were reported 
for the 1.025 in. tube with the veiocity at disappearance 
of Regime III being about 3.5 ft./sec. In the case of the 
1.5 in. tube Regimes I, II and III are clearly apparent up 
to water velocity of 2.50 ft./sec., although only at that 
water velocity were sufficiently high air rates realized to 
illustrate Regime IV. Insufficiently high air rates were 
obtained at the water rate of 3.47 ft./sec. definitely to 
determine the presence of Regimes II, III or IV, although 
the correlation to be described later suggests that Regime 
III is actually missing and that Regimes II and IV are 
merged. With the 2.50 in. tube Regimes I, II and III are 
apparent to water velocities of 2.05 ft./sec., although 
Regime III becomes increasingly indistinct as the water 
velocity is increased. It is clear that the disappearance 
of Regime III and the mergence of Regimes II and IV 
takes place at different water and air velocities in the 
tubes of different diameter. 

The general shape of the holdup curve is more or 
less the same for all diameters and all superficial water 
velocities. The holdup ratio is unity (slip velocity is 
zero) at zero air-water ratio, it increases with air-water 
volume ratio passing through an inflection point about in 
the middle of Regime I, and thereafter increases to a 
second inflection point near the end of Regime II, after 
which it increases more slowly, levels or decreases. 


°A more detailed form of this paper has been deposited as Document 
No. 5723 with the ADI Auxiliary Publications Project, Photoduplication 
Service, Library of Congress, Washington 25, D.C. A copy may be 
secured by citing the Document No. and by remitting $1.25 for photo- 
prints, or $1.25 for 35 mm. microfilm. Advance payment is required. 
Make checks or money orders payable to:' Chief, Photoduplication Service, 
Library of Congress. 
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the Test Section Length 22.88 Ft. 
: feet Test Section Mid Point Pressure 36.0 Psia 
cu Test Section Mid Point Temperature as Indicated 
: Specific Volume of Liquid Phase = 0.01603 Ft’/Lb 
nn 5 Specific Volume of Saturated Gas Phase as Indicated 
hese Superficial Water V eloci ities Approximately Constant at 0.850—0.873 ft/sec. 
rates == == | = = SSS i => 
‘hich (1) (2) (3) (4) (5) (6) (7) | (8) | (@) (10) (11) | (12) 
S the —_-- SS ee es F _ aa — a | ri i * a, } ‘ = 
lumn | . ; er be es Specific 
Tube Water | Saturated | Mid Point | Discharge Test | Volume Pressure | Pressure 
> test Test | Diameter Rate Air Rate |Temperature} Air-Water | Section | Saturated | Holdup Drop Drop Flow 
gives Number | inches ft®/sec ft? sec* °F Ratio** | Air-Water | Gas Phase |Ratio***|ft H.O/ft | Regime | Pattern 
test Ratio | ft/lb 
nap, 280D 0.630 | 0.00183 | .00426) 62 2.33 1.29 | 5.37 1.81 | 0.413 | I Slug 
i 281D 0.630 | 0.00183 | .00574 | 62 3.14 1.73 5.37 1.82 | 0.348 | I Slug 
6 to 282D 0.630 | 0.00183 | .00828 | 62 4.52 2.12 5.37 2.13 | 0.300 | I Slug 
» Bes 283D 0.630 | 0.00183 loo | 63 8.47 2.85 5.38 StS 0.300 II Slug 
. 284D 0.630 | 0.00183 .0230 | 63 | 12.6 2.79 5.38 4.52 0.330 II Froth 
irmng 85D 0.630 | 0.00183 0281 | 64 w4 | 2n | 32 5.68 0.360 II Froth 
of 86D 0.630 | 0.00183 | .0334 | 64 18.2 | 3.02 | 5.39 6.03 0.372 | III Froth 
and 87D 0.630 | 0.00183 | .0407 | 64 22.2 | 3.33 5.39 6.66 | 0.376 | III Ripple 
288D 0.630 0.00183 .0308 64 16.8 2c 5.39 6.20 0.372 II Froth 
0; 289D 0.630 0.00183 .0451 64 | 24.6 4.49 | 5.39 5.48 0.388 IV Film 
_ 290D 0.630 | 0.00183 | .0576 | 64 31.5 4.91 5.39 6.42 | 0.425 | IV Film 
“ize 
rate 380D 1.025 | 0.0050 | 0.0620 | 69 12.4 2.32 5.45 5.35 | 0.312 | II Froth 
dro 381D 1.025 0.0050 0.117 | 69 23.4 3.65 5.45 6.41 0.288 | III Froth 
P 382D 1.025 | 0.0050 | 0.158 70 31.6 4.34 5.47 7.28 | 0.293 | IV Ripsie 
(b) 383D 1.025 | 0.0050 | 0.137 75 | 27.4 | 4.13 5.53 6.63 | 0.295 | III Ripple 
vater 386D 1.025 0.0050 0.0303 73 6.06 2.47 $.3% | 2.45 0.256 | I Slug 
drop 387D 1.025 0.0050 0.0568 | 74 1 1 4 2.26 S251 | §.04 0.330 | II Froth 
arly 388D 1.025 | 0.0050 0.0775 | 72 Ioe2 2.81 5.51 §.5§2 0.321 III Froth 
a 380D 1.025 | 0.0050 | 0.0157 | 73 3.20 1.83 5.51 1.75 | 0.346 | I Slug 
€ in 390D 1.025 0.0050 | 0.0285 | 72 5.80 2.28 dca 2.54 0.272 I Slug 
7 to 391D 1.025 0.0050 | 0.0383 73 7.60 2.74 5.51 2.04 0.252 I Slug 
the 392D 1.025 0.0050 0.00788 | 70 1.60 0.930 5.47 1.72 0.498 I Slug 
» | = ~ c ce ec . 
> Ail 345S 1.50 0.0107 0.396 79 St .2 5.02 Saas 7.41 0.217 III Ripple 
.630 346S 1:50 | 0.0107 | 0.214 | 66 20.1 3.52 5.41 5.71 | 0.281 III Froth 
ater 350S 1.50 0.0107 | 0.282 | 67 26.5 3.97 5.42 6.68 0.244 Ill Froth 
rible 3518 1.50 0.0107 0.170 67 16.0 2.70 5.42 5.93 0.320 III Froth 
352S 1.50 0.0107 | 0.128 67 12.0 2.37 5.42 5.06 0.350 III | Froth 
ared 353S 1.50 | 0.0107 0.0898 67 8.43 2.46 5.42 3.42 0.337 | II Froth 
the 3548S 1.50 0.0107 0.0321 | 67 3.01 1.47 5.42 2.05 0.377 I Slug 
rted 3558 1.50 0.0107 | 0.0587 | 67 Dads 2.04 5.42 2.70 0.320 | I | Slug 
nce 356S 1.50 0.0107 | 0.0197 | 79 1.84 0.914 5.55 2.01 0.540 | I | Slug 
. 2578 1.50 | 0.0107 | 0.0493 | 79 4.63 1.81 5.55 2.56 | 0.293 | I | Slug 
the 258S 1.50 | 0.0107 | 0.0888 | 80 8.34 2.70 5. 56 3.09 | 0.335 | II | Froth 
¢ up 259S 1.50 0.0107 0.338 81 31.7 4.94 3.57 6.42 0.241 III | Ripple 
that 360S 1.50 |0.0107 | 0.389 80 36.5 veh 5.56 7.09 0.232 Ill | Ripple 
I to 481S 2.50 | 0.0296 | 0.0586 | 78 1.98 0.900 5.54 2.20 | 0.433 | I Slug 
vere 4825 2.50 0.0296 | 0.0731 | 79 2.48 0.992 $55 2.50 0.481 I Slug 
’ to 483S 2.50 | 0.0296 0.0940 | 79 3.20 1.12 aio 2.85 0.432 I Slug 
ugh 484S | 2.50 | 0.0296 | 0.126 79 4.25 1:25 yoo ise 0.377 I Slug 
2 285S 2.50 | 0.0296 | 0.163 78 5.50 1.41 5.54 3.90 0.340 I Slug 
pine 286S | 2.50 | 0.0296 0.199 78 6.74 1.53 5.54 4.41 0.330 I Slug 
are 487S 2.50 | 0.0296 0.250 78 8.45 1.74 5.54 4.86 0.340 II Froth 
are 488S 2.50 | 0.0296 | 0.336 78 11.3 1.98 5.54 5.70 | 0.297 III Froth 
uch 489S 2.530 0.0296 0.295 79 10.0 1.89 S59 5.30 0.328 Ill Frot 
8 490S | 2.50 0.0296 0.423 79 14.3 2.22 S59 6.45 0.277 Ill Froth 
ater 491S | 2.50 0.0296 0.510 78 17.2 2.47 5.54 6.95 0.284 III Ripple 
nce _ | 
IV — - 
the * At test section temperature and average pressure. 
** Column 4 divided by column 3. 
or **Column 6 divided by column 7. 
ater 
Ries Interpretation and correlation of results locus exists for the first minimum point or that the air- 
ater Effect of diameter on flow pattern water volume ratio at which the first pressure drop mini- 
r in ‘ As has previously been shown the flow patterns may mum and the transition from Regime I to Regime II 
dD 8 best be defined and correlated by reference to the pressure occurs is independent of tube diameter. In the cases of 
fter drop regimes which in turn are defined by the two minima the maximum point and the second minimum point, how- 
and the maximum points in the unit pressure drop-air- ever, the loci depend upon diameter. The locus of the 
water volume ratio curve. The loci of the abscissa values maximum points and that of the second minimum points 
tion of these points as functions of the superficial water intersect at an increasing water rate and a decreasing air- 
ae ’ velocity are shown in Figure 19. This figure presents the water volume ratio as the tube diameter increases. In- 
ired. data for all tube diameters including the data for the 1.025 sufficient data were obtained to completely define the 
vice. 8 : os e I a 
: in. tube reported in Part I. It will be seen that a common maximum and second minimum loci for the 1.50 and 2.50 
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in. tubes but the trends are quite well established and are 
indicated by dashed lines. 

These loci lines are further correlated and related to 
the flow patterns in Figures 20, 21, and 22. Figure 20 
shows the single locus of the first minimum points for all 
four diameters and the flow-pattern transition lines. The 
bubble-slug transition appears as a single line for all 
diameters. The slug-froth transition lines are slightly 
different for each diameter because of the effect of dia- 
meter on the locus of the maximum points. The different 
loci for the maximum points were, by trial, found to be 
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brought closer together by plotting superficial water 
velocity V,, vs. RyD% as in Figure 21. The spreading of 
these loci at high water velocities reflects the effect of 
water velocity on their mergence with the loci of the 
second minimum points. On this plot the slug-froth transi- 
tion again appears as a family of lines, one for each 
diameter. Figure 22 presents the locus of the second 
minimum points as functions of V;, vs. RyD-%. Again the 
abscissa was determined by trial to bring the three separate 
loci lines of Figure 19 into coincidence. The transitions 
from froth to ripple and from ripple to film flow each 
appear pretty well as single lines on this plot. 
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Effect of diameter on holdup 

The holdup data of Figures 2-18 inclusive and that 
presented in Part I were cross plotted and analyzed to 
determine the effect of tube diameter on holdup. By trial 
it was found that all the holdup data for Regimes I and Il 
could be correlated in terms of the superficial water 
velocity and the parameter RyD™%. The correlation is 
indicated in Figure 23. The data for Regime III correlate 
well in terms of the parameter Ry alone as shown in 
Figure 24. The data for Regime IV were found to corre- 
late reasonably well with the parameter RyD~% as illus- 
trated in Figure 25. 








Effect of diameter on pressure drop 

Che total pressure drop was separated into its two 
1/3 . . ~ ~~ * . 
ie components by calculation using Eq. (1). Friction factors 
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were then evaluated from Eq. (2). Numerical values of 
the pressure drop components and the friction factors are 
presented in Table 2 deposited as Document No. 5274* 
with the American Documentation Institute Auxiliary 
Publication Project, Photo Duplication Service, Library 
of Congress, Washington 25, D.C. The general nature 
of the effect of tubing diameter on the friction factor is 
shown in Figure 26 which presents data for all four 
diameters at a typical constant water velocity of some 


°A more detailed form of this paper has been deposited as Document 
No. 5724 with the ADI Auxiliary Publications Project, Photoduplication 
Service, Library of Congress, Washington 25, D.C. A copy may be 
secured by citing the Document No. and by remitting $1.25 for photo- 
prints, or $1.25 for 35 mm. microfilm. Advance payment is required. 
Make checks or money orders payable to: Chief, Photoduplication Service, 
Library of Congress. 


201 





















’ 
L 





FRICTION FACTOR, f 





— emcrteeccinaren — on 
© —V,=0.869 ft. per sec., D*2.50 in. 
X -V,=0.869 ft. per sec., D*1.50in. 
& —V20.873ft. per sec., D=1.025in. 


© —V,*0.852ft. per sec. 0=0.630in. 
a | | 


| 


2.0 50 20 50 joo 200 
AIR-WATER VOLUME RATIO, Ry ,cu. ft. per cu. ft. 


Figure 26—Typical friction factor relationships. 


0.87 ft./sec. By trial it was found that a very good corre- 
lation was obtained by plotting friction factor vs. DV, 

with a parameter of V,,D*°, where Vq is the superficial 
gas phase velocity bend upon the average flowing tem- 
perature and pressure and upon the total tube cross section, 
ft./sec. This is shown in Figure 27. The abscissa DV,, was 
chosen as representing a reduced form of the superficial 
Reynolds number for the liquid phase. The parameter 

V, D2. ® may be considered the product of a reduced form 
iw the superficial Reynolds number based on the gas phase 
and the term D2-5. The conventional single phase friction 
factor corresponding to a parameter value of zero, is 
shown for reference purposes. The two phase friction 
factor curves lie well above and to the right of the single 
phase line but approach the latter at high abscissa values. 


Conclusions 


The new data and correlations enable the prediction, 
for any diameter between, say 0.5 and 3.0 in., of the flow 
pattern, the holdup and the pressure drop as functions of 
the diameter and the air and water rates for the upwards 
vertical flow of air-water mixtures at an average pressure 
and temperature of 36.0 psia and 70-85°F. respectively. 
The correlations are restricted in application to systems 
where the average pressure is constant at 36.0 psia and 
where the fractional change in specific volume of the 
water saturated gas phase is small enough to justify the 
use of an average valve. The effect of average pressure 
(or gas phase density) is now under study and will be 
reported upon shortly. 
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A Heat Exchanger Design 


for Heavy Water Reactor Service’ 


A. D. DUFF, JR.2 


The volume of heavy water required for a power 
reactor plant using this moderator has an important 
influence on investment. Likewise, the cost of leakage 
justifies extraordinary containment measures. A heat 
exchanger design is described which satisfies these 
requirements, as well as those imposed by the radia- 
tion field in which the equipment must operate. 


EAVY water heat exchangers in the primary coolant 
H loops of a nuclear power plant pose many design 
problems not encountered in conventional practice. These 
problems arise from three basic considerations as shown 
in Table 1. First is the high cost of heavy water, currently 
quoted by the United States Atomic Energy Commission 
at $28.00/lb. Second is the radioactivity induced in the 
heavy water stream. Third is the necessity for maintain- 
ing the heav y water at a very high level of i isotopic purity 
in order to maintain its efficacy as a moderator in the 
reactor. 

TABLE 1 


| NUCLEAR HEAVY WATER 


Basic DESIGN CONSIDERATIONS | 
| HEAT EXCHANGER 

@ Cost of Heavy Water — $28/Ib. 

@ Radioactivity in Heavy Water 

@ High Isotopic Purity Requirements 


Owing to the cost of heavy water, the hold-up in the 
exchangers represents a substantial investment and 
requires that optimization calculations include the cost, 
normally neglected for conventional fluids. This mater- 
ially affects the usual ratio of volume to effective heat 
transfer surface. Relatively high over-all heat transfer 
coefficients and friction losses are a natural consequence. 
The volume contained in connections from the piping 
system to the heat exchangers again has an important 
effect on design of heads and connecting piping. The 
cost of heavy water also requires a high degree of integrity 
— leakage in the design. For example, a leak rate of 
5 drops/min. represents an annual loss of about $10,000.00. 

Radioactivity in the heavy water arises from tritium 
formed by exposure of deuterium to neutrons, from 
activity induced by neutron bombardment of corrosion 
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products and other trace impurities and from fission 
products in the event of fuel failure. This condition lends 
added weight to the need for freedom from leakage, 
because of the potential hazard to personnel resulting 
from entry of radioactive heavy water into the secondary 
loop, for ‘example. Similarly, maintenance of radioactive 
equipment is both time-consuming and expensive, because 
of limits.on exposure of personnel, the need for protective 
clothing, and the like. The heat exchanger, therefore, must 
be so designed as to permit reduction of radioactive levels, 
by cleaning, to a point where reasonable maintenance work 
can be done, at intervals, over a relatively long period. 
This requires that it be possible to decontaminate the equip- 
ment without disassembling it and, preferably, without 
necessitating close contact between maintenance personnel 
and the equipment. Ideally, the exchanger should be free 
from places where contaminants can lodge, and should be 
capable of being decontaminated by some relatively 
simple technique, such as a high volume flush, with 
provision, of course, for disposal of the contaminants. 
Minimizing the necessity for performing any main- 
tenance work at all is fully as important as simplifying 
maintenance functions. This is, of course, true in con- 
ventional plant design, but, for the reasons stated earlier, 
assumes far greater proportions for the nuclear case. 


The need for maintenance of isotopic purity of the 
moderator is still another reason for designing against 
leakage. The entry of secondary coolants, such as light 
water, into the heavy water stream requires installation 
and operation of separating equipment, which can repre- 
sent an appreciable cost. 


The factors mentioned so far tend to make the cost 
of heat transfer surface for the nuclear heavy water 
exchanger higher than for normal utility or chemical 
_ service. The achievement of tolerable costs for 
large units makes it important to use available materials 
in sizes and specifications reasonably close te normal 
industrial practice. Similarly, consideration must be given 
to feasibility of fabrication of exchangers in large units 
with minimum necessity for developing techniques and 
skills which depart w idely from those employed in the 
fabricating industry. 


As summarized in Table 2, the heavy water heat 
exchanger in a nuclear plant primary coolant loop must 
be designed for: 

(1) Low heavy water hold-up. 

(2) Minimum leakage in either direction. 

(3) Minimum maintenance. 
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(4) Ease of required maintenance. 
(5) Minimum fabrication cost considering the above 
requirements. 


TABLE 2 

NUCLEAR HEAVY WATER 
DESIGN REQUIREMENTS 

HEAT EXCHANGER 
Low Heavy Water Holdup 
Minimum Leakage in Either Direction 
Minimum Maintenance 
Ease of Required Maintenance 
Minimum Fabrication Cost 


The remainder of this paper will describe how an 
actual design was developed and executed to meet these 
requirements, together with comments on fabrication 
experience and on the ultimate service performance of 
the equipment. The facilities described were designed for 
transferring heat from a large flow of heavy water to a 
light water stream, but the principles involved are appli- 
cable to other secondary media. 

The parameters of heat load, temperature, and pres- 
sure were first determined, after which several alternate 
basic designs of oe were considered. One concept 
(shown in Figure 1) w yhich may be called the plate type, 
consisted of sheets 7 either stainless steel or aluminum, 
pressed into a multi-channel form. A stack of such sheets 
would be assembled, and welded or brazed at points of 
contact. Manifolding arrangements at the ends would 
permit flow of light water and heavy water through 
alternate channels. This concept had ‘the advantage of 
efficient use of material, good heat transfer coefficients 
on both sides, low cost, low heavy water volume, and 
probably fairly good decontaminability. Some small ex- 
perimental units were fabricated and ‘tested, but it was 
found that necessary fabrication techniques would require 
considerable dev elopment. Furthermore, leak repair, either 
during fabrication or as a result of failure in service, 
appeared impractical, and it was judged that the degree 
of approach to perfection in fabrication technology 
required to attain zero leakage was beyond any reasonable 
expectation. 

Heat exchangers based on the falling film principle, 
using either horizontal or vertical tubes, were next con- 
sidered. Calculations and tests showed, however, that the 
hold-up for this design was surprisingly high, by com- 
parison with other types. An important factor in holdup 
for this design was the requirement for dissipating the 
high inlet velocity which resulted from the use of com- 
paratively small pipe in the main loop. 


A number of variations of the shell and tube type of 


2 


Figure 1—Plate type heat exchanger. 


exchanger were then considered. The requirement of low 

ratio of volume to heat transfer surface favored the 
smallest practical tube diameter, with heavy water inside 
the tubes. This led to consideration of 4 in. and 4 in O.D. 
tubes fastened into tube sheets by rolling, welding, 
drifting, or combinations of these. Inv estigation rev ealed 
that it was phy sically possible to expand or weld tubes 
of these sizes into tube sheet holes. However, the problem 
of controlling such operations to produce the required 
tightness of joints, considering the thin tube wall and the 
delicacy of tools small enough to enter the tubes, made 
the use of such small tubes unattractive. 


An alternate approach (indicated in Figure 2) to the 
use of these very small tube sizes was suggested by Alco 
Products Company and consisted of a number of tube 
ends welded or brazed together into a pack, with a ring 
welded to the periphery of the pack and a reducing 
fitting to the ring. The reducer, in turn, would be fastened 
to a tube sheet by conventional means. Diameters of 1 in. 
to 2 in. were visualized for that end of the reducer which 
would be fastened to the tube sheet. It was visualized 
that the tubes after leaving the end cluster would be 
splayed out and held in a ‘rectangular pattern by small 
square tube supports. Experimental assemblies of this type 
were made and leak-tested. Considerable difficulty was 
experienced in securing “Freon”-tight joints. The diffi- 
culty lay in securing a deep enough joint between 
adjacent tubes to make an effective seal. A further draw- 
back was that there were no means for correcting leaks 
after assembly of tubes into the pack. Thus, there would 
be a substantial amount of heat transfer surface lost if 
leaks showed up after fabrication, since it would be 
necessary to plug the entire pack. 

Since no practical method of using small tubes had 
been found, it was concluded that the best approach lay 
in the use of larger tubes individually fastened to tube 
sheets, with heavy water volume minimized by placing 
cores in the tubes. Discussions of the problems were held 
with manufacturers, who, at first, were inclined to regard 
2 in. O.D. tubes as the smallest practical size for this 
equipment. Numerous calculations and some development 
work led to the use of 3} in O.D., type 304 stainless steel 
tubes, with 0.3 in. diameter solid core rods. Double tube 
sheets were considered essential to permit rapid detection 
and location of any leakage which might occur, to mini- 
mize loss of heavy water into the light water stream, and 
to collect leakage, thereby permitting longer operation 
under leakage conditions. Use of maximum length tubes 


SECTION "A-A coven 34 


Figure 2—Tube Pack heat exchanger. 
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was settled upon so as to minimize number of units and 
joinrs, thereby reducing cost as well as number of 
potential sources of maintenance. 

One interesting aspect of the heat transfer calcula- 
tions may be briefly mentioned. Tube side flow was 
found to lie in the transition region between laminar and 
well-defined turbulent flow. Because of the unusual ratio 
of length to diameter, the calculated heat transfer co- 
efficient was found to be somewhat lower than might be 
expected. Both tube side and shell side coefficient were 
verified experimentally to minimize allowances for un- 
certainty. 

For the final design, tube size was selected as } in. 
O.D., 18 gage, type 304L stainless steel. Tubes were to 
be 29 ft. 6 in. long which was the greatest length available 
from most producers. Design features of the unit, such as 
number of tubes, core size, bafHing, and the like, were 
established by optimizing calculations described in a paper 
by Brinn and Cichelli™. 

In carrying out the tube sheet design, it was calculated 
that a thickness of 7 in. to 9 in. would be required if 
unstayed sheets were used, because of the large shell 
diameter. It was soon apparent that the cost and time 
required for drilling these tube sheets would be intoler- 
able, even if it were assumed that the problem of con- 
trolling drift of the drill over this length could be solved. 
The tube sheets were therefore recalculated using stay- 
bolts and a thickness of 1} in. found adequate. In com- 
puting stresses, credit was taken for support of the outer 
tube sheet by the inner one via the tubes, since the 
staybolts did not bear on the inner tube sheets, as shown 
in Figure 5. 

The shell (Figure 5) was of conventional cylindrical 
fixed tube sheet design. Stainless steel plate was used, 
rather than carbon steel, so as to avoid expansion joints 
which were believed to be a potential source of main- 
tenance. The shell closure was effected by welding it 
to the inner tube sheets. Tube supports, cross flow baffles, 
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Figure 3—Heat exchanger in place with the inlet nozzle 
covered for testing. 
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impingement baffles, and inlet and outlet nozzles were 
conventional. 

The heads were gasketed to the other tube sheet, as 
shown in Figure 5, by means of rectangular cross section, 
neoprene gaskets, which were essentially fully confined 
and compressed by controlled bolting. Consideration was 
given to using a welded joint, but it was decided that 
conservative design and meticulous fabrication could 
produce an essentially maintenance-free gasketed joint. 
Stay-bolts were threaded into the outer tube sheet and 
passed through bosses on the heads with nuts on their outer 
ends which served to transfer the tube sheet load to the 
head. O-rings on the bodies of the bolts provided a seal. 
As mentioned previously, the exchanger shell terminated 
at the inner tube sheet at each end with the space between 
the two tube sheets, in each case, being closed by a 
removable shroud. ; 

The head design, shown in Figure 4, is an interesting 
feature of this heat exchanger and was described in a 
paper by Cichelli & Boucher (2). The basic requirements 
were low volume, low pressure drop, and reasonably equal 
distribution of fluid to the tubes, in addition to furnishing 
structural support to the tube sheets. The problem was 
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Figure 4—Heat exchanger head details. 
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Figure 5—Heavy water heat exchanger. 
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complicated by the numerous staybolts used for support- 
ing the outer tube sheets, which lie across the water flow 
path. Figure 3 shows the head in place with the inlet 
nozzle covered for testing. 

Each heat exchanger was supported on a pair of stand- 
ard railroad car trucks. These facilitated handling within 
the building and guided the unit as it expanded and 
contracted in service. In practice, one set of trucks was 
anchored and the other allowed to float. 


Mention was made earlier of the desirability of util- 
izing established techniques and standard materials. It is 
worth noting that even though this principle was followed, 
a substantial portion of the engineering effort required for 
these heat exchangers was devoted to fabrication and 
material procurement. Perhaps the largest procurement 
effort was that of obtaining the necessary tubing in the 
required quantities. The efforts of several redraw mills 
were needed to produce the required amounts, and, of 
course, it was necessary to standardize testing and dimen- 
sional control so that tubes from the several mills could 
be used interchangeably. Seamless, butt-welded, and so- 
called swaged tubing all were used, and experience was 
good with all types. 


An essential feature of the fabrication phase was ex- 
haustive testing. After each row of tubes was inserted 
and fastened, tube to tube sheet joints were leak-tested, 
first with “Freon” and then with hydrostatic pressure. It 
was found important to make the “Freon” test prior to 
the hydrostatic test, because it had been demonstrated 
earlier that water could effectively seal apertures through 
which the “Freon” would otherwise indicate a leak. The 
first exchanger was also thermally cycled to simulate 
operation and uncover any foreseen effects of temperature 
change. 

As a step to minimize eventual contamination of the 
system in which the heat exchangers were to operate, the 


units were cleaned by shooting wool felt pellets throu; gh 
the tubes with a compressed air gun. Heads were cleaned 
separately by conventional methods. 


Actual experience with operating these exchangers has 
been very favorable. Heat transfer coefficients have been 
found close to predicted values. A few leaks have oc- 
curred in service at head gaskets, and staybolts, and in 
two cases, tubes chafed through on tube supports. How- 
ever, the number of such faults has been extremely low 
by any standard. The extreme measures taken to assure 
freedom from leakage or from the need for maintenance 
have been profitable. Build-up of activity in these ex- 
changers has occurred, however, with the result that 
some recent process changes made were relatively costly 
due to the limitations on exposure of personnel. In design- 
ing new exchangers for similar service, even greater atten- 
tion to details for preventing entrapment of radioactive 
particles would seem warranted. Development of tech- 
niques for the use of small tubes could result in substantial 
cost reductions. 


In conclusion, the authors feel that the initial require- 
ments of low volume, minimum leakage, minimum main- 
tenance, and ability to decontaminate were achieved at a 
cost which was reasonable in the light of the requirements. 
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First Order Rate Processes and Axial 
Dispersion in Packed Bed Reactors 


JAMES J. CARBERRY? 


By employing recently secured data for the axial 
dispersion of mass flowing through fixed beds, it is 
shown that for first order isothermal processes con- 
ducted in fixed beds of practical aspect ratios (L/D,), 
the effect of mixing upon kinetics is insignificant. It 
follows that in shallow bed systems usually employed 
in kinetic research efforts, the influence of axial 
mixing cannot be ignored, particularly in the case 
of a liquid reactant. The dispersion data and the 
analyses support the contention of Ergun “'*) and the 
recent treatment of the problem by Epstein “°). 


9 importance of fixed bed processes in the chemical 
and allied industries has motivated numerous theo- 
retical and exper rimental analyses of the distribution of 
heat and mass in fixed beds. Early attention was given to 
the problem of radial diffusion of mass by Bernard and 
Wilhelm |), In a recent text, Smith (2) discusses the work 
of those who concerned themselves with the radial mixing 
problem. 

The attention focused upon the distribution of heat 
and mass is obviously justified by the simple fact that 
point reaction rates in a packed bed depend uniquely 
upon concentrations of heat and mass at a given point. 
Obviously, the goal of rational reactor design demands a 
basis w hereby these concentration distributions may be 
predicted. The ingenious design method of Baron (?) 
illustrates this point emphatically. 


Axial dispersion 


The problem of radial dispersion having been treated 
rather thoroughly, that of axial or longitudinal dispersion 
has only recently attracted attention. Danckwerts ‘), 
Kramers, and Alberda “), and Carberry and Bretton ‘®) 
studicd the dispersion of tagged water flowing through 
fixed beds. McHenry and Wilhelm ‘7) determined axial 
diffusivities for gas systems while Carberry and Bretton 
(6) also studied the axial dispersion of a gas in the Reynolds 
number region below that considered by McHenry and 
Wilhelm “), 

The general method of correlation is that of relating 
the Peclet number, D,v/E, to Reynolds number, D,v p/ tt. 
The data of those inv estigations cited above are shown in 
Figure 1, based upon average bed velocity. Discussions of 
these studies are given elsewhere ‘7: §), 

The limiting v value of the Peclet number for axial 
dispersion in fixed beds has been theoretically shown to 
be about two by Aris and Amundson °®), McHenry and 
Wilhelm (7) and Carberry “'). If a fixed bed is viewed 


1Manuscript received May 16, 1958. 
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as a series of mixers, then, as stated by Kramers and 
Alberda (5) as well as others, the numbers of mixers n is 
related to the bed length L, fluid velocity v and dispersion 
coefficient E, by 

me WIRE cde 


When perfect mixing prevails in each void cell of 
length yD,, where y is the ratio of void cell length to 
particle diameter (<1), then n = L/yD, and 


1/yD, = v/2E or 2/y = D,v/E = Peclet........(2) 


This is a condition achieved at high levels of turbulence, 
corresponding to a void cell mixing efficiency of unity. 
A model relating cell mixing efficiency to the Peclet 
number has been presented earlier (1°). The model shows 
that cell mixing efficiency, ¢ is equal to (y/2)Pe. 

The limiting value of the Peclet number has been 
predicted and verified experimentally for a gas system 
by McHenry and Wilhelm ‘) while the date of Carberry 
and Bretton ‘®) for water appear to approach this limit at 
high Reynolds numbers as shown in Figure 1. 


Axial dispersion and reaction 


Hulburt “), Danckwerts (4), Wehner and Wilhelm (12) 
and Yagi ‘!8) have treated mathematically the case of 
first-order kinetics and axial dispersion. While Wehner 
and Wilhelm “'2) have discussed the subtleties of the 
boundary conditions, it is agreed that the fraction of 
material unreacted, f, at the bed exit is related to the dis- 
persion and chemical kinetic variables by the equation ‘*) 


f = 4a/[(1 + a)? exp (— (Pe’/2) (1 — a)) — 
(1 — a)? exp (— (Pe’/2) (1 + a)].......(3a) 


where a = Vi + 4k0/Pe’.... jason 


k@ = first-order time constant 
Pe’ = Peclet number in terms of bed length, Lv/E 


He-AIR $ 
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Figure 1—Peclet number vs Reynolds number. 
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Since n = Lv/2E, then Eq. (3) may be written 


f, = 4a/[(1 + a)? exp(— n(1 —a)) — 
(1 — a)? exp (— n(1 + a)) 


Vi + 2k0/n.. 


When piston or plug flow prevails, n +> and 


where a 


fon = exp (— kO)....... tS) 
while as n + 1, 
f,; = 1/(1 + WO isssx ; ; (6) 


The ratio f,/f,, thus represents the fraction of material 
unreacted in the effluent of a fixed bed isothermal reactor 
equivalent dispersionwise to n inixers as compared to the 
fraction unreacted for a bed of infinite mixers (E — 0). 
This ratio is plotted in Figure 2 as a function of k@ or 
kL/v. Various levels of conversion for a piston flow 
reactor are also noted. 


Recalling that since 


n = Lv/2E = (D,v/E) (L/2D, RG 2s 5 cna 





ke 


Figure 2—Relative influence of axial mixing upon first- 
order processes. 


and noting the range of Pe values found experimentally 
and shown in Figure 1, it is immediately apparent that 
axial dispersion exerts only a trivial effect upon a first 
order kinetic reaction conducted isothermally in a fixed 
bed of practical dimensions. This is particularly evident 
for gaseous systems as the McHenry and Wilhelm data 
(7) demonstrate that virtually perfect void cell mixing 
(Pe = 2) is achieved beyond ; a Reynolds number value 
of dent 30. While the effect ha liquid systems is more 
pronounced, it is evident that a low \ alue of Pe is not 
critical so long as the aspect ratio, L/D,, remains high. 
Thus the dispersion data and kinetic analysis simply 
rationalize the common observation that only shallow 
beds exhibit back-mixing characteristics. 

With respect to shallow bed systems, the work of 
Ergun ‘'#) deserves notice. Analyzing systems employed 
by laboratory research investigators, i.e., shallow beds 
and/or differential reactors, Ergun demonstrated that in 
mass transfer studies, piston flow does not alw ays prevail 
and in consequence, the use of the log mean driving force 
to evaluate transfer coefficients is inadmissible. The reader 
is referred to Ergun’s papers ‘'4.'5) for details. Since dis- 
persion data were not available at the time of his analyses, 
Ergun could only consider the extremes of piston flow 
and perfect mixing throughout a bed. His point would 
seem to be deserving of notice insofar as the evidence 


now at hand would suggest that deviations from piston 
flow certainly may be critical i 


experiments carried on 
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in small scale fixed bed reactors where low aspect ratios 
(L/D,) are not uncommonly encountered. It is clearly 
of some significance for an investigator to know w hether 
a given research reactor behaves as 2, 10 or 50 mixers, 
since rate coefficients can only be accurately derived from 
a true knowledge of the reactant concentration profile 
between bed inlet and exit. For ex cample, Epstein “'%) did 
not ignore this point in analyzing mass transfer data 
secured in shallow beds. He was able to demonstrate, 
using the idealized model employed by McHenry and 
Wilhelm (7?) to explain their gas dispersion data, that the 
assumption of piston flow inherent in the usual integration 
of the mass transfer equation led to transfer coefficients 
which were only slightly in error for the gas system under 
study. 

Epstein has developed a correction factor, the deriva- 
tion of which is to be published by him, which permits a 
modification of the log mean driving force in terms of 
the number of perfect 1 mixing lengths, L/D,, — 
ing dispersion behavior when Pe = 2 and + y = 1. This 
relation is of obvious generality since for any ae of 
the Peclet number the number of equivalent mixers n 
equals (L/D, )e or Lv/2E. Defining R as the ratio of inlet 
to outlet drivi ing forces, Epstein’s correction factor is ‘'®), 


Y = InR/n(R!" — 1). (8) 


The use of this relationship with the data of Figure 1 
should prove to be a useful guide to those concerned 
with the procurement of accurate rate coefficients in 
laboratory scale fixed bed systems. 

For the gas mass transfer system considered by Epstein, 
it was shown, for R = 8 and n = 8, that Y equalled 0.88. 
This value of n corresponds to a Pe value of 2 in Epstein’s 
gas system. For a liquid system at the same Reynolds 
number, the value of Pe would be about 0.5; equivalent 
to n = 2. With the same driving force ratio, R, the 
correction factor becomes 0.57. Thus the value of the 
transfer coefficient differs by 75% from that calculated 
by neglecting dispersion in a liquid mass transfer system 
for the same driving force ratio. 


The contentions of Ergun ‘!5) on mixing in fixed beds 
appear worthy of notice “with respect to ‘liquid systems 
while Epstein’s analysis “'®) of gas-solid networks at 
Reynolds numbers greater than 30 is supported by the 
dispersion data. 


\t extremely low values of the Reynolds number for 
<1) the data of Carberry and Bretton ‘*) shown 
in Figure 1 suggest that the axial dispersion coefficient is 
indeed equal to the molecular diffusivity of the gas 
system. In consequence, the Peclet number decreases 
linearly with decrease in Reynolds number in this region 
of flow. Research studies conducted therein will certainly 
demand consideration of the dispersion effects upon the 
rate process under study. 


gases ( 


No doubt one could re-examine the research literature 
and find numerous instances where a serious mixing 
problem has been ignored in evaluating rate constants 
and/or } factors. It suffices to cite the work of Gamson, 
Chodos and Hougen ‘'7) who studied the rate of evapo- 
ration of water from pellets in an air stream. Bed depths 
not oe 24 in. were used for particle sizes up to 
} in, It is obvious from this discussion that at best, only 
aa three mixing stages prevailed in the study with the 
largest particle size. Resnick and White ('8) studied mass 
transfer in fixed beds of } in. depth. Aspect ratios of about 
12 were encountered in their system, however, very low 
values of the Reynolds number were —— Since the 
systems studied comprised naphthalene in air, CO, and H, 
respectively, the actual Peclet numbers iamecunaiian the 
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systems may well have been in the transition region 
= een perfect cell mixing and that dictated by mole- 

lar diffusivity. Indeed Resnick and W hite (18) report 
an inexplicable effect of particle diameter. As we have 
seen, for a given bed depth, particle diameter determines 
the dispersion behavior of the system. 


\ final point of some importance is that of entrance 
effects in shallow beds. The author ‘8) found that for a 
given Reynolds number the dispersion coefficient, FE, in- 
creased significantly as bed length is reduced to 6 in., an 
event anticipated by Ergun“) some time ago. This 
anomaly is attributed to acceleration effects as the fluid 
passes from an unpacked section into the bed. These data 
merely emphasize the importance of situating dummy beds 
immediately upstream of an active bed in any shallow 
bed research work. 


Conclusions 


In the light of recently correlated axial dispersion 
data for fixed beds, it is evident that for first order 
isothermal processes carried out in packed beds of prac- 
tical aspect (L/D,) ratios, the effect of axial mixing on 
the kinetics is negligible. 

On the other hand, in bench scale studies in fixed 
bed reactors, often characterized by low aspect ratios, 
caution must be exercised in the evaluation of transfer 
or rate coefficients since the assumption of piston flow 
may not be justified for liquid systems. As shown by 
Epstein (1), the mass transfer equation based upon piston 
flow behavior can be applied, even at low aspect ratios, 
to a gas system where each void cell acts as a perfect 
mixer (Pe = 2), provided that a suitable correction 
factor is applied to the log mean driving force. At low 
values of the Reynolds number shallow bed experiments 
in gas or liquid systems demand a consideration of axial 


mixing as it affects concentration gradients within the 
system. 
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Nomenclature 


a = function defined by equations 3b and 4b 
D, particle diameter, cm. 
S axial dispersion coefficient, cm.?/sec. 
fraction of reactant remaining at bed exit 
first-order rate constant, sec.~ 
natural log 
bed length, cm. 
number of perfect mixing lengths, defined by equation 1 
Peclet number in terms of bed length 
= Peclet number in terms of particle diameter 
= ratio of inlet to outlet driving force as defined by Epstein 
16) in equation 8 
average fluid velocity in fixed bed 
log mean driving force correction factor ‘'®), equation 8 
exponential, e 
= nominal holding time, L/v, sec 
ratio of void cell length to partic le diameter 
void cell mixing efficiency 


Subscripts 


one mixer 
n mixers 
infinite number of mixers 
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Correction Factor for Axial Mixing 
in Packed Beds’ 


NORMAN EPSTEIN? 


Based on the model of a fixed bed as a series of 
perfect mixers, a factor is derived for correcting the 
logarithmic mean driving potentials in the conven- 
tional plug-flow equations for fluid-particle heat and 
mass transfer in packed beds. This factor, which 
accounts for axial mixing and is a measure of the 
transfer efficiency of the bed relative to the idealized 
case where axial mixing is absent, increases towards 
unity with bed height but decreases with ratio of 
terminal driving potentials. 


ECENT experimental studies 1; 2: 3. 4) on axial mixing in 
fixed beds have been rationalized by viewing a fixed 

bed as a number of perfect mixers in series ( % 5°), 
According to this model, the number of perfect mixers 
in a fixed bed 1s given by 
LV 


2E AI) 


The axial Peclet number Pe has experimentally () and 
theoretically ‘4. 5.6) been found to reach a limiting value 
of about two for turbulent gas flow through a bed of 
spheres. For other cases, values smaller than two have 
been obtained (: * The number of perfect mixers in 
a packed bed is thus equal to or less than the number of 
packing layers L/D,. 

For either heat transfer or mass transfer between a 
fluid stream and a packing surface maintained at uniform 
temperature and/or concentration, the assumption of no 
longitudinal mixing leads to the conventional integrated 
rate equations containing the logarithmic mean of terminal 
driving potentials. The actual presence of axial mixing in 
a packed bed precludes the use of these plug- flow equa- 
tions, unless a suitable correction factor F is employ ed to 
modify the log mean driving potential. Such a factor 
would be analagous to those used to modify LMTD in 
multi-pass shell-and-tube heat exchangers ‘, The modi- 
fied equations would then appear as 


q = hA (At)wF (2) 
for heat transfer, and 
N = kA (AY)imF (3) 


for mass transfer. The derivation of an expression for F 
follows. 


Derivation 


Consider a fluid entering a packed bed at temperature 
t,, transferring heat to the particle surface which is main- 
tained at uniform temperature t*, and leaving the bed at 


iManuscript received August 7, 1958 
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temperature t,. Such would be the case, for example, in a 
fixed bed of particles undergoing gas drying within the 
constant rate period (8). If the bed is viewed as a series 
of n perfect mixing cells each having surface area A/n 
and constant heat transfer coefficient h, then for the first 
mixing cell, steady state heat transfer results in 
WCp (t. — ti) = h (A/n) (t: — t*)... 
Rearranging Eq. (4), 
t. + Pt* 


where 


~ WC,n 
Similarly for the second cell, 
ti + Pt* 
1+ P 
Combining Eq. (6) and (4a), 


to + 2Pt* + P?*t* 
(1 + P)? 


t= 


E = (1 + P)? 
to .° 
Similarly for the third cell, 

+ Pt* 


ts = 
1+P 


Combining Eq. (7) and (8) 


+ 3Pt* + 3P%* + P*t* 
ts 
(1+ P)3 


se 1+ Pp) 
= (1 - ) 
ts ¢” 


This pattern continues to the nth mixer, resulting in 


t. — t* 
= (1+ P).. . (10) 


t* 


that is 
R = (1+ P)*.. ; «ist ROD 


P= Rm —] (10b) 
\n overall heat balance on the packed bed yields 


q = WC, (to — ta).. (11) 
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first 


(4) 


(4a) 


Comb ‘ning Eq. (11), (2) and (5), 
to — tn 
(At )im 7 


» : . . ‘ * 
But for uniform surface temperature t*, 


o — t*) — (ta — t* 
(At)im 2 (t ) ( Es _ 
t. — t* 


tn — t* 
Therefore Eq. (12) is equivalent to 
InR 


Combining Eq. (14) and (10b), 
In R 
n (R!/2 — 1) 


Eq. (15) can be used in conjunction with Eq. (2) and 
(11) to predict integral heat transfer rates in a packed 
bed. 
The derivation of F for mass transfer follows exactly 
—- lines to that for heat transfer and again vields 
. (15), where R in this case is defined as the ratio 
a * — Y,)/(Y* — Y,). The mass balance analagous to 
Eq. (11) is 
N = G(Y, — Y,).. (16) 


Eq. (15) and (16) can be used in conjunction with Eq. 
(3) to predict integral mass transfer rates in a packed bed. 


Discussion 


The axial mixing factor F may be considered a measure 
of the heat or mass transfer efficiency of a packed bed, 
being unity in the absence of back-mixing. It has already 














F= AXIAL MIXING FACTOR 
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been applied by the author in a previous paper ‘*). Though 
rigorously Eq. (15) applies only where t* or Y* is uni- 
form throughout the packing, it will probably give a 
reasonably good answer in Eq. (2) or (3) for many cases 
of non-uniform surface conditions, especially if a judi- 
cious average value is used to represent the surface 
condition. 

Figure 1 is a plot of Eq. (15). It is seen that where the 
number of perfect mixers is small, that is, for short packed 
beds, axial mixing markedly reduces the transfer efficiency 
of the bed. This is particularly true if the longitudinal 
temperature or concentration change of the fluid is 
large, resulting in a high value of driving potential ratio 
R. In addition, end siteane in shallow beds probably 
increase the axial diffusivity, thereby decreasing n accord- 
ing to Fq. (1), and hence the efficiency of transfer is 
further decreased. The limiting case is that where n = 1, 
in which the entire bed acts as a perfect mixer. 

On the other hand, for deep beds, where n is large, 
axial mixing can usually be neglected in determining 
integral transfer rates. The limiting case here is an 
infinitely long bed, for which F = 1.0. 

Mitigating against both extremities is the fact that, in 
practice, a shallow bed is usually accompanied by a low 
potential drop and therefore low value of R, while a deep 
bed is usually accompanied by a high value of R. 

Due to neglect of F in Eq. (2) and (3), many reported 

values in the literature of particle-fluid sera average heat 
and mass transfer coefficients for packed beds are prob- 
ably low, especially where shallow beds have been used. 
Discrepancies berween results of various investigations are 
probably in many cases also caused by this neglect. An 
attempt to reconcile some of these discrepancies by 
means of Eq. (1), (2), (3) and (15) is currently being 
undertaken by the author. This is being supplemented by 


experiments on very shallow beds. 


Nomenclature 


= surface area of packing, sq. ft 


heat capacity of fluid at constant pressure, B.t.u./(Ib.) 


(°F.) 


= particle diameter, ft 


axial diffusivity, sq. ft./hr. 
axial mixing factor, dimensionless 
flow rate of non-diffusing fluid, lb.-moles/hr. 
= particle-fluid heat transfer coefficient, B.t.u./(hr.) (sq. ft.) 
(°F.) 
particle-fluid mass transfer coefficient, lb.-moles diffusing 
component/(hr.) (sq. ft.) (AY) 
= depth of packing, ft. 
number of perfect mixing cells in series, dimensionless 
mass transfer rate of diffusing component, lb.-moles/hr, 
hA/WC,n for heat transfer, dimensionless 
ky A/Gn for mass transfer, dimensionless 
axial Peclet number = D,V/E, dimensionless 
heat transfer rate, B.t.u./hr. 
ratio of inlet to outlet driving potential: (t, — t*) 
(tn — t*) for heat transfer, (Y* — Y.,)/(Y* — Yj) for mass 
transfer, dimensionless 
temperature of fluid stream, °F. 
surface temperature of packing, °F. 
t — t", F. 
interstitial fluid velocity, ft./hr. 
weight rate of fluid flow, lb./hr. 
molal ratio of fluid stream, Ib.-moles diffusing component/ 
Ib'-mole non-diffusing fluid 
molal ratio at packing surface, lb.-moles diffusing com- 
ponent/Ib.-mole 
AY Y* — Y, Ib. moles diffusing component/Ib.-mole non- 
diffusing fluid 


Subscripts 


OQ = packed bed entrance 
Bey; n Ist, 2nd, 3rd, nth perfect mixer 
Im = logarithmic mean of terminal conditions 
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Some Procedures for the Evaluation 
of Reactor Fuels and Sheathing Materials 
at Chalk River’ 


R. F. S. ROBERTSON? 


In a power reactor the fuel will be operating 
under very severe conditions. Consequently before a 
fuel element can be designed it must be tested under 
conditions which are similar to those in the reactor, 
and the effect of a large number of variables on the 
behavior of the fuel must be investigated. 

This paper discusses the fuel development pro- 
gram at Chalk River with special reference to the 
testing of uranium oxide reactor fuels in loops 
inserted in the NRX reactor. In these loops water 
may be recirculated over the test elements at tem- 
peratures up to 550°F. and pressures up to 2000 
psi. The characteristics of a loop are described and 
the results of various tests are briefly discussed. 

The testing of sheathing materials such as Zir- 
caloy-2 and aluminum-nickel alloys is also described. 


HE main objective, towards which most work at Chalk 

River is directed, is the production of economic 
nuclear power. There are many routes to this objective, 
some apparently more promising than others, and the one 
which Atomic Energy of Canada Limited has chosen is 
the use of natural uranium dioxide as the fuel and heavy 
water as moderator and coolant. 

This paper does not pretend to be in any way a 
presentation of technical data, but rather is a general 
discussion of the program to examine the irradiation 
behavior of reactor fuels and sheathing materials which is 
now underway at Chalk River. At the end of the paper 
will be found a list of recent AECL publications which 
give more technical details. 

Figure 1 shows a sch<matic representation of the 
proposed reactor NPD-2 which is being built in con- 
junction with the Canadian General Electric Co. Limited 
and Hydro Electric Power Company of Ontario. It can 
be seen that in essence there are two separate heavy water 
circuits—a high and a low temperature circuit. The high 
temperature and pressure circuit is one in w hich the cool- 
ant (D,O) is recirculated past the fuel which is in the form 
of UO, rods sheathed in Zircaloy-2, a corrosion resistant 
alloy of zirconium and tin. The fuel is held in roughly 
130 channels, each channel consisting of a tube made of 
Zircaloy-2, which is capable of constraining the high 
pressures at the temperatures involved. The heat which 
the heavy water picks up from the fuel is removed in the 
heat exchangers where water (H,O) is boiled and the 
resultant steam is passed to the turbines. The cooled D,O 
returns from the heat exchanger to the pumps and thence 
back to the reactor. 


1Munuscript received July 25, 1958. 

2Senior Research Officer, Chemistry and Metallurgy Division, Atomic 
Energy of Canada Limited. Chalk River, Ont.—presently on loan to 
Argonne National Laboratories, Lemont, Ill. 

This article is based on a paper presented at the Joint A.I.Ch.E.-C.L.C. 
Chemical Engineering Conference, Montreal, April 20-23, 1958. 
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Figure 1—Schematic representation of proposed NPD-2 
power reactor. 


The second, low temperature and pressure circuit is 
known as the moderator system. The D,O is held in a 
large aluminum vessel which surrounds the pressure tubes. 
It is hag rong at a relatively low temperature (about 
120°F.) by circulation through heat exchangers. This 
sy a although it has its own special problems, will not 
be discussed in this paper which will deal only with the 
problems associated with the high temperature circuit. 
The discussion will be further limited to the problems 
associated with the fuel and sheath. Allison“) discusses 
some of the problems associated with the interaction of 
the high temperature coolant with the rest of the system. 

Table 1 gives some information on the NPD-2 reactor, 
the fuel and some of the conditions under which it will 
be operating. 


TABLE 1 


MATERIALS AND OPERATING CONDITIONS FOR A POWER REACTOR 


Materials 

Fuel Material UO, 

Density 10 g/cc (92°; theoretical) 
Sheath Zircaloy-2 
Zircaloy-2 


Pressure Tubes 


Stainless steel, carbon steel, 
Inconel 


Other materials in high 
temp. circuit 


Operating Conditions 

80,000 kw 

20,000 kw 
1,000 psi 


Thermal Power 

Electrical Power 

Pressure 

D,.O Temperatures 
Entering Reactor i. 
Leaving Reactor F. 


Heat flux (fuel surface to water) 135,000 (BTU) (hr!) (ft-?) 


Thermal Neutron Flux 5 x 10'3 (n) (em?) (see!) 
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It is quite apparent that the fuel in NPD-2 and in any 
similar power reactor, will be operating under quite 
stringent conditions, and moreover under conditions re- 
garding which there is but little experience. Thus it would 
be a waste of time to design, and unwise to build such a 
reactor before information was available regarding the 
behavior of the fuel elements under conditions as nearly 
as possible those expectcd in the reactor. 

To this end an extensive testing program is proceeding 
at Chalk River and experiments are being conducted both 
the reactor and in the reactor. This paper will 
discuss only in-reactor experiments. 


out of 


Experimental loops 


known as “Loops” have been 
inserted in the NRX reactor at Chalk River. Westing- 
house Atomic Power Division, Pittsburgh, U.S.A. have 
been responsible for the insertion of the original high 
temperature and pressure loops in NRX and these loops 
have been operated in close cooperation with them. A 
schematic representation of a typical loop is shown in 
Figures 2 and 3. It can be seen that a loop i is really a small 
version of the primary coolant circuit of a power reactor 
and in some ways can be looked on as a pilot plant. 


Several test facilities, 


DELAYEO wWEUTROW 
wOomiTOR 


COOLER 


r i | 


10m 
EXCHANGE 
BY-PASS 
SEE 
Fig.3 


HEATER 


TEST 
SECTION 


2 — Schematic diagram of Loop X-2 in the NRX 


reactor. 


Figure 


Figure 3—lon exchange by-pass in Loop X-2. 


One of NRX 
reactor and in Its place a stainless steel tube is inserted. 
The 


this 


the uranium rods is removed from the 
1 and removed from 
The coolant, 
the 


the canned 


inserted 
the 


test can be 
the 


fuel unde 


tube through closure at 


top. 
water in bottom 
and led out 
roto! 
of 


sources of heat, the heater and the heat produced i 


most cases, 1s brought in through 


near the top. It is recirculated by 


and is heated electrically to 


the 


pump a temperature 
S00-S50°F. by heater. Thus there are two possible 


n the 
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The 


test specimens themselves. Temperatures are maintained 
at a fixed level, whether the reactor is operating or not, by 
a judicious juggling of the power supplied to the heaters 
and the flow rate of coolant to the coolers. 

To prevent local boiling of the coolant, the loop is 
pressurized by the surge tank which is merely a large 
vessel containing water and steam heated to a suitable 
temperature to provide the required pressure. The surge 
tank also acts as a reservoir to accommodate changes in 
the volume of the coolant caused by temperature changes, 
In addition, by initiating a flow of coolant through the 
vapor space of the surge tank and thence back to the loop, 
the surge tank may be operated as a very efficient degasi- 
fier. 

On a bypass are found the ion-exchange resins which 
are used to maintain the purity of the water. Because 
these resins decompose at high temperatures the loop 
water must be cooled before it reaches them. 

Table 2 shows typical conditions which can be main- 
tained in a loop. 


TABLE 2 


Loor OPERATING CONDITIONS 


Pressure Up to 2,000 psi 


Coolant Temp. Up to 570°F. 


Linear Flow Past Fuel Specimens —_ Up to 20 (ft.) (sec.~') 


Maximum thermal neutron flux ~ 4x'10'(n) (em>?) (seem!) 


A variety of sizes and shapes of fuel may be irradiated 
in the test section. Rods up to 8 ft. long by 3 in. in 
diameter may be irradiated or the same length of fuel 
can be subdivided into a “sausage chain” of short rods 
joined end to end. Flat plates have also been irradiated. 
Most work to date has been done on short rods of up 
to a foot long and roughly one half inch in diameter. 

At present the UO, fuel is made by pressing UO, 
powder into cylindric al pellets at high pressures. These 
pellets are then sintered in a hydrogen atmosphere for 
several hours at a temperature about 1700°C. The 
resultant ceramic has a density of about 10.5 gm. (i.e. 
95° theoretical). Fuel elements are made by 
1e proper size 


of 
cc. 


about OF 


slipping the pellets into Zircaloy tubes of t' 
and welding the ends closed. 


Fuel examination facilities 


Having irradiated the under the correct 
conditions it is of course, imperative that they be examined 
to see how they have behaved. Thus the Remote Control 


Handling F acility, or “C: is a Most important comple- 


specimens 


ave” 
ment to the loop program, 

It is possible to remove samples from the loops and 
them cave 
room 


take shielded containers, to the metallurgy 
which is shown in Figures 4 and 5. This is really a 
surrounded with thick concrete walls into which are set 
thick lead glass windows. Remote handling manipulators 
pass through the walls and thus the specimens which are 
highly can examined Table 3 
shows some of the equipment which is used in the cave. 


radioactive be remotely. 


PABLE 3 


\PPARATUS IN CAvi 


Stereo microscope and camera, (magnification 1x to 90x) 
\icro-manipulator 

Vec hanic al testers | Hardness, 
Protilometer 

Vex hank al h it k 
Lathe 

Micro drill press 


Grinding and polishing wheels (for me tallographic work ) 


Send Strength, Impact, Tensile) 


SAW 
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Figure 4—Photograph of remote handling cell. 
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Figure 5—Vertical section of remote handling cell. 


Experimental investigations 

The discussion of experiments carried out at Chalk 
River using these facilities will be limited to the testing 
of UO, and sheathing materials suitable for use in future 
water cooled power reactors. 


1. Irradiation behavior of UO, 


Of course one of the most important questions is: 
“How does the oxide stand up to irradiation?” In this 
connection it must be remembered that because a large 
amount of energy is being released in the fuel by the 
fission process and because dense U QO, is a relatively poor 
conductor of heat, very large temperature gradients will 
occur throughout the fuel. The centre may well be at a 
temperature rather close to the melting point of 2800°C., 
whereas the surface of the element will be less than 500°C, 
This is a temperature gr: adient of about 9000°C. per inch. 
Many specimens of UO, have been irradiated in order to 
determine the effect of high thermal neutron flux and 
the consequent high temperature gradient, on the fuel 

In most instances UO, has withstood the irradiation 
very well, and beyond some cracking, caused probably by 
thermal shock, the pellets appear visually and dimension 
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Figure 6—Two samples of compacted and sintered UO:. 

Both irradiated to about 4800 MWD/Tonne UO:.. One has 

cracked and the other is uncracked. (1 Tonne = 10° 
kilograms) 


Figure 7—Compacted and sintered UO... irradiated to 
about 4800 MWD/Tonne UO. shows formation of central 
flue. 


ally unchanged even after quite long irradiations. Figure 
6 shows a typical example. 

In some samples, in which the core temperature has 
been over an estimated 1600°C., grain growth at the centre 
of the pellet has been seen. In other samples, of composi- 
tion UO,,, the interesting phenomena shown in Figure 
7 has been seen where a large central “flue” has been 
formed down the centre of the specimen, and long 
crystals radiate from this flue to the edge of the specimen. 
There is still much speculation as to the exact nature of 
the process occurring in these specimens, some believe it 
to be melting while others believe it to be a mass transport 
phenomenon. 

Howey er, it is now possible to tell the reactor designers 
that as long as certain conditions are met (and some of 
these conditions can be outlined quite specifically), UO, 
is a fuel which will withstand long irradiations w ithout 
swelling or distorting in any way. 


2. Fission product gases 


Among the products of fission are stable isotopes of 
the rare gases Kr and Xe. At the high temperatures found 
in the fuel, atoms of these gases can diffuse from the 
oxide lattice and collect in the space between the fuel and 
the sheath. After a very long irradiation it is possible that 
enough of these gases may have collected in this space 
to give high pressures inside the sheath and consequent 
deformation or rupture. To study this question, specimens 
of fuel in Zircaloy sheaths are irradiated and after irradi- 
ation the sheath is punctured in a specially constructed 
apparatus in the cave and the amount of gas released is 
measured, 

To date this problem of pressure build-up does not 


appear to be a serious one from a practical standpoint, 
but much information is being obtained regarding the 
actual mechanisms of release and subsequent history of 


these vases. 


1958 215 





3. Defect tests 
Information is being obtained on the release of fission 
product gases and also of other fission products by irra- 
diating samples in which the Zircaloy sheath has been 
punctured by drilling a small hole through it to the fuel 
below. By analysing the loop water, it can be determined 
which isotopes escape most readily, the extent to which 
they escape and also their subsequent history in the loop. 
In other words, do they remain in solution, do they end 
up on the ion exch: ange resins or do they deposit on the 
walls of the loop piping? 
Studies such as these give information on several other 
problems such as: 
(i) Is the oxide stable to high temperature water and to 
steams 
(ii) What will the problems be in an operating reactor 
if the fuel sheath ruptures and what necessary steps 
must be taken for decontamination? 
(iii) What is the most suitable monitoring system to give 
early warning of sheath failure? 


1. Corrosion studies 

To date an alloy of zirconium and tin, known as 
Zircaloy-2 appears to be the most promising sheath for 
enclosing UO, fuel because it combines desirable nuclear 
properties Ww ith good resistance to corrosion by high 
temperature water. Figure 8 shows the exterior of a 
Zircaloy sheated specimen after an exposure in loop 
water at a temperature of 500°F. for several months. The 
original machining marks are clearly visible, and obviously 
no corrosion has taken place and no fouling of the surfaces 
has occurred. 


: 


zt 


Figure 8—End of Zirealoy-2 sheathed fuel specimen — 
exposed to water at 525°F. for four months. 


‘There are probl ms associated with Zircaloy, however. 
One problem is that, as well as using Zircaloy as a sheath 
material for the fuel, it is proposed to use it in the pressure 
tubes of NPD-2, which must have a lifetime of at least 
tn years. When Zircaloy is in contact with water at 
500°F. a dark highly protective film of ZrO, forms over 
the surface resulting in a low corrosion rate ( 1.4 x 10°55 
in/yr) which is independent of such factors as water 
velocity etc. However after a certain time which is not 
definitely known but will probably be about eight years 
under NPD water conditions, the protective oxide layer 
changes to a white less protective oxide. The corrosion 


rate is still low (~2.5 x 10° in./yr.) but whether it is 
now dependent on such factors as water velocity is not 
known for certain. Tests are now being conducted in 
which the transition of black to white oxide has been 
induced to occur in a short time by treatment of the 
metal in steam at 750°F. These specimens, with the transi- 
tion film, are now being tested with water at 500°F, 
moving past them at 20 ft./sec. to see whether the resultant 
corrosion rates are acceptable. 

Although Zircaloy-2 is a good sheathing material, it is 
expensive. The search for a cheaper sheathing material 
is continuing and a number of alloys of aluminum with 
nickel are being studied. These alloys have quite low 
corrosion rates (~0.003 in./yr.) in water at 500°F. pro- 
vided certain conditions are fulfilled. This rate is just 
acceptable for a fuel sheath provided that the alloy is 
not susceptible to pitting and providing it is not necessary 
to irradiate the fuel for too long a time. The processes 
occurring are not yet fully understood but it is hoped 
that by the addition of corrosion inhibitors to the water 
corrosion rates can be decreased. 

In this paper only the highlights of the AECL testing 
program have been considered. Appended is a list of 
pertinent AECL reports, all of which are unclassified, 
which give fuller details of various aspects of the program. 


Reference 


(1) Allison, G. M., C.J.Ch.E. 36: : 220: st. 58. (AECL Report 
DCI-30). 


AECL Reports on Loop Tests, etc. 


Booth, A. H., “A Method of Calculating Fission Gas Diffusion from vo, 
Fuel and its application to the X-2-f Loop Test.”” AECL Report 
CRDC-721 (AECL No. 496), Sept. 1957. 

Booth, A. H., Bain, A., Rymer, G. T., Kerr, D. B., “Collection and 
Measurement of Stable Fission Xenon in Sheath Puncture Experiments 

Method and Results to Date.” AECL Report CRDC-719, Oct. 1957. 

Booth, A. H., “A Suggested Method for Calculating the Diffusion of 
Radioactive Rare Gas Fission Products from UO, Fuel Elements.” 
AECL Report DCI-27, Sept. 1957. 

Tomlinson, M., Robertson, R. F. S., Pattison, J., ““The Chemistry of the 
first AECL Fuel Rod Test in the WAPD-CR-IV-X2 Loop.” AECL 
Report CRDC-650, June 1957. 

Robertson, R. F. S., Krenz, F. H., “Defect Test on U-2 w/o Zr Alloy 
in the X-2 Loop. Test No. 1.”” AECL Report CRDC -646 (AECL No. 
366), April 1956. 

Robertson, R. F. S., “Defect Test on U-2 w/o Zr Alloy in the X-2 Loop. 
lest No. 2.”” AECL Report CRDC-647. 

Robertson, R. F. S., G. M. Allison, “Irradiation of Compacted and 
Sintered UO, in the X-2 Loop Tests X-2-f and X-2-i.”” AECL Report 
CRDC-718, May 1958. 

Robertson, R. F. S., “Review of Loop Programme in the NRX 
Reactor.”” AECL Report UK/C-6 LOL, Sept. 1957. 

Howieson, J., “A Report on the First NPD Single Element Irradiation 
and Failure.”” AECL Report CRNE-744 (AECL No. 552), August 
1957. 

Bain, A. S., Payne, W. E., Howieson, J., “Irradiation of UO, and U-2 
w/o Zr Alloy in the CR-V Loop.’”” AECL Report NEI 84, Feb. 1958. 

Bain, A. S., “Construction and Shielding Properties of Remote Control 
Handling Cell, Metallurgy Building, Chalk River, Ontario.” AECL 
Report CRL-39, March 1957. 

Krenz, F. H., “The Corrosion of Aluminum-Nickel Type Alloys in High 
femperature Aqueous Service.”” AECL Report CRL-40, May 1957. 

Krenz, F. H. et al, “Chalk River Experience with Zircaloy-2 and Alumi- 
num-Nickel-Iron Alloys in High Temperature Water.”” Paper No. 
P/194 submitted to the Second United Nations International Con- 
ference on the Peaceful Uses of Atomic Energy, Geneva, Sept. 1958. 

Chalder, G. H. et al, “Fabrication and Properties of Uranium Dioxide 
Fuel.” Paper No. P/192 submitted to the Second United Nations 
International Conference on the Peaceful Uses of Atomic Energy, 
Geneva, Sept. 1958. 

Robertson, J. A. L. et al, “Behaviour of Uranium Oxide as a Reactor 
Fuel.” Paper No. P 193 submitted to the Second United Nations 
International Conference on the Peaceful Uses of Atomic Energy, 
Geneva, Sept. 1958. 
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Conditioning of D,O 


in Heavy Water Power Reactors 


G. M. ALLISON? 


This report indicates the problems associated with 
the maintenance of the purity of the D,O in a 
pressurized heavy water power reactor and how these 
problems are handled. 

In a pressurized water power reactor, impurities 
will be introduced into the coolant and moderator 
from the corrosion of the steel system. These impur- 
ities are mainly particulate and colloidal in nature. 
They must be kept to a minimum since when circu- 
lated through the reactor they become activated by 
the neutron flux and when deposited throughout the 
system they make accessibility for maintenance diffi- 
cult. Impurities may also be introduced into the water 
from fuel element ruptures. These also create main- 
tenance problems since they are intensely radioactive. 


Heavy water systems require a constant guard 
against the introduction of light water which would 
degrade the D,O and make the operation of the 
reactor inefficient or perhaps impossible. The isotopic 
content of the heavy water must be checked cons- 
tantly. Upgrading may be required periodically or 
continuously. 

The general purity of the water in the reactor is 
maintained by passing a fraction of the main circu- 
lating volume through a by-pass containing a filter 
and a mixed bed ion exchange resin column. This 
column serves to remove ionic impurities as well as 
acting as an efficient filter itself. It is also used to 
maintain the pH of the water at the desired level 
whether neutral or alkaline. In the case of heavy 
water systems the resin must be deuterized before use 
to prevent degrading the D,O. 


A PRESSURIZED heavy water power reactor essentially 
consists of a closed system in which heavy water 
coolant is heated by the fuel rods and is circulated 
through heat exchangers and back to the reactor vessel 
again. Boiling is prevented by pressurizing the system. A 
simple diagram of such a system is shown in Figure 1 of 
reference 1. The out-reactor parts of this system are 
mild steel and the in-reactor pressure tubes are made of 
Zircaloy-2 (a zirconium-tin alloy ). 

There are several ways in which impurities can be 
introduced into the D,O. 

(1) The main impurities in the system will be the 
corrosion products (so-called “crud”) which are formed 
in the reaction between the hot water and the metal 
surfaces. This crud will be present mainly in colloidal 
and particulate forms and to a very slight extent as ionic 
mat.rial. These are the same corrosion products that one 
finds in a conventional steam plant. 


tManuscript received July 25, 1958. 

2Associate Research Officer, Chemistry and Metallurgy Dvision, Atomic 
Energy of Canada Limited, Chalk River, Ont 

This article is based on a paper presented at the Joint A.LL.Ch.E.-C.1.C 
Chemical Engineering Conference, Montreal, April 20-23, 1958. 


Fuel elements showing deposited crud after irradiation for 
one month in pH7 water at 525°F., conductivity < 1 umho 
em'. 


(d) 


(c) Fuel element after irradiation for one month in pH7 
water, 2 x 10-* M in K.SO,, at 525°F. 
(d) Fuel element after irradiation for two months in 
pH10 water at 525°F. 


Figure 1 


(2) Fission products will be present from the fissioning 
of uranium occurring as an impurity in the Zircaloy-2 
and as so-called “tramp” uranium on the surface of the 
sheathed fuel. This tramp uranium is very difficult to 
remove completely after fabrication of the fuel elements. 
Fission products will also be introduced into the coolant 
when fuel element failures occur. 


(3) Light water if introduced into the system will 
degrade the D,O. 
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The crud problem 


In any power reactor employi ing light or heavy water 
at high temperatures corrosion products will cause certain 
problems which are the same as those in a conventional 
steam plant and others which are characteristic only of 
nuclear plants. In steel systems, crud is formed by cor- 
rosion reactions and is composed mainly of magnetite, 
Fe,O,, with small amounts of the oxides of the alloying 
metals in the steel, Cr, Ni, Mn, etc. The number and 
amount of the alloying elements appearing in the crud 
will depend on the composition of the steel, mild steel 
crud being nearly pure Fe,O,. If oxygen is present in the 
water Fe,O, will also be present. 

Some of the corrosion products will remain on the 
metal surface as a corrosion film while some will be 
releas:d to the water and carried throughout the system 
mainly in particulate and colloidal form. In a conventional 
steam plant this material can be detrimental in many ways. 

(a) The particulate crud can build up in places causing 
reduction in flow rate. 

(b) It can cause wear and malfunction of moving 
parts made to close tolerances. 

(c) It can deposit on heat transfer surfaces. 

Nuclear plants also have these problems and in addition 
they have two more— 

(d) Under certain water conditions, crud will prefer- 
entially deposit on the fuel elements causing reduction in 
heat transfer rates and possible over- heating of the fuel 
with subsequent burnout. 

(e) Corrosion products become activated by absorp- 
tion of neutrons in the reactor and are deposited through- 
out the out-reactor parts of the system. These activation 
products are y-emitters and many have sufficiently long 
half-lives to cause maintenance problems. 

Pressurized water nuclear power plants then will have 
more problems due to corrosion products than conven- 
tional steam plants. It must be remembered that in nuclear 
plants amounts of crud which would be considered 
insignificant in a conventional steam plant can cause 
serious problems due to radioactivity and to preferential 
deposition on the fuel elements. 
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Figure 2—Resin deuterization curve. 


It is thus desirable to keep the crud concentration in 
the water as low as possible and to choose water conditions 
which will also prevent preferential deposition on the fuel, 
Considerable data on the behavior of crud in stainless and 
mild steel systems containing high temperature water at 
high pressure have been obtained from the loops described 
by Robertson “) where simple flow diagrams of a typical 
loop are shown in Figures 2 and 3. It has been shown quite 
definitely in loop experiments in the NRX reactor at Chalk 
River and in American reactors, that maintaining the high 
temperature water at pH 9.5 - 10.5 with KOH or LiOH 
keeps the crud level in the water low, prevents the prefer- 
ential deposition of the corrosion products on the fuel 
elements and removes any crud which may have been de- 
posited under former, less desirable, water conditions (2: 8), 
It has also been shown that neutral water which has been 
made 2 X 10-* M in potassium sulphate has the same effect 
as high pH water “). The effect of potassium sulphate has 
been shown in one loop test only while that of KOH has 
been substantiated many times. More K,SO, experiments 
are required. Table 1 contains results obtained in loop 
experiments in NRX in which different water conditions 
prevailed. The effect of high pH and K,SO, on the crud 


level in the water can be readily seen. 


TABLE 1 


EFFECT OF WATER CONDITIONS ON CruUD LEVELS IN Loops 


Water Condition Crud Level, 


ppb* 


pH 7, conductivity < 1 puke (Sealalees Steel) “10- 50 
pH 7, 10-4 M K2SO, (Stainless Steel ) 1-5 
pH 10, KOH (Stainless Steel) 2-10 
pH 10, KOH (Mild Steel) 10—50 


*ppb = parts per billion. 


Large variations in the crud level are often caused by 
thermal cycling of the loops from reactor trips and start- 
ups and from deliberate cooling of the loops for main- 
tenance purposes. Temperature changes tend to crack 
the corrosion film causing crud to be released to the loop 
water. In one of the mild steel loops the crud concentra- 
tion under high pH water conditions has risen to over 
1000 ppb for short periods due to thermal cycling. 

Photographs of Zircaloy-2 clad UO, fuel elements 
which were irradiated under conditions given in Table 1 
are shown in Figure 1°), The effect of water conditions 
on the deposition of crud on the fuel can be readily seen. 


Radiation problems 


(a) Active crud: Besides the physical problems per- 
taining to crud which have just been considered, there is 
also an activity problem. Corrosion products when ex- 
posed to a neutron flux in the reactor vessel become 
radioactive. The extent of the activity will depend on the 
absorption cross section of the material for neutrons, the 
length of time in the neutron field and the half-life of 
the radionuclide formed. 

Active corrosion products found on pipe walls out- 
side the reactor vessel can come from three sources: 

(i) from corrosion of active metal in the reactor vessel, 
the active crud being carried throughout the system in 
the circulating water, 

(ii) from activation of crud formed outside the reactor 
but deposited within the reactor vessel for a time and 
then transferred throughout the system, 

(iii) from crud being activated as it is repeatedly 
carried in and out of the reactor in the recirculating 
water. 
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Tne last source is insignificant compared to the other 
two. 

Table 2 shows the major corrosion product activities 
foun! in crud samples collected from the recirculating 
water in a stainless steel loop in NRX. 


TABLE 2 


Y-EMITTING ACTIVATED CORROSION PRODUCTS 


% of total y- 

Nuclide| Half- Mode of Formation activity in crud 
Life sample after 7 

days 


45d | neutron capture by Fe*® 5-15 


| 

| Only a major 

| contributor for the 
hy first 24 hours 


2.6 hr | neutron capture by Mn® 


7id neutron capture by Ni®§* 
5.2 yr | neutron capture by Co®? 
27d | neutron capture by Cr*® 

2.8d | neutron capture by Sb’?! 


60d | neutron capture by Sb!** | J 





(Zr%® + Nb®) 





65d | neutron capture by Zr*! | 40-50 


*In this case neutron capture is followed by a proton emission 
instead of a y-ray. 


Although antimony is a very minor constituent in 
steels it has a large absorption cross-section for thermal 
neutrons and antimony activities have been found to be 
major contributors to the long-lived y-activity in crud 
collected from the water in both stainless and mild steel 
loops in NRX. They can also be formed by a (p,n) 
reaction on the tin in the Zircaloy. Zr®® comes from the 
activation of zirconium in the Zircaloy-2 used to clad the 
fuel elements and used as structural material for the pres- 
sure tube in this particular loop. Nb® is the daughter of 
Zr®® and is also a y-emitter (half-life 35 d). 

It can readily be seen that a build-up of long-lived 
y-emitting active corrosion products on pipe walls could 
cause maintenance problems due to parts of the system 
becoming inaccessible except for short periods. The maxi- 
mum permissible total body irradiation which anyone is 
allowed is 5 roentgens (r)* per year which is roughly an 
average of 100 milliroentgens (mr) per week. A man 
working in a general radiation field of 1 r/hour therefore, 
will receive his week’s dosage in 6 minutes. As we will 
see later, high fields are experienced in loops when "G 
element ruptures occur. From the activity stand-point, i 
is thus important to keep the crud release rate and a 
corrosion rate low by using proper water conditions to 
minimize the transfer of crud throughout the system. 


(b) Activation products of water: The water in a 
pressurized water reactor will become very radioactive 
due to neutron activation of O' to form N' (half-life 
7.3 seconds), and O' to form N' (half-life 4.1 seconds). 
If the water is not deaerated, natural A*°® will be present 
and be activated to form A‘! (half-life 1.8 hours). Some 
F's (half-life 1.9 hours) is also formed from O78 by a 
(proton, neutron) reaction. Fortunately, these activation 
products are short-lived and do not cause radiation prob- 
lems shortly after the reactor shuts down. They do, how- 
ever, necessitate shielding of piping containing cooling 
water circulating between the reactor and heat exchangers. 


®°The roentgen is a frequently used unit of radiation dose. 


(c) Fission product activities: Fission products will 
be present in the cooling water in a power reactor due 
to uranium impurities in the fuel sheathing material and 
pressure tubes and due to “tramp” uranium on the sheaths. 
However, the major source of fission products will be 
from ruptured fuel elements, i.e. where the sheaths fail. 


Loops in NRX have been operated with deliberately 
defected fuel elements in order to study the release of 
fission product activity from different types of fuel and 
to study the behavior “of the different fission products in 
the cooling water. By passing a portion of the main loop 
flow through a mixed bed ion exchange column, a much 
lower equilibrium activity level in the loop water can be 
attained for the longer- -lived fission products (ti, > a 
few hours) than if no ion exchange is used at all. If the 
fission product source is removed, the water can be 
cleaned up satisfactorily by the resin column. 

Some fission product activities tend to deposit on pipe 
walls with the corrosion products. These are mainly the 
more insoluble elements. Table 3 shows the major y-emit- 
ting active nuclides found deposited on pipe walls in a 
loop in NRX when defected UO, elements were present. 
In this loop the in-reactor pressure tube was stainless steel. 


TABLE 3 


Major y-EMITTING NUCLIDES DEPOSITED ON PIPE WALLS IN 
DeFEctT TEST 


Relative % of 
Total y after 
4 days 


Source 


Niclide Half-Life 


S. Crud 
Prod. 


Mo x 


Sb!24 


Co” 
Fe59 
Ru 103 
Tel29 
Tels 
Ba!0 
Cr! 
Zr95 


The piece of pipe removed from the loop for analysis 
was Type 304 Stainless Steel, 1 in. Schedule 80, and had 
had loop water flowing through it at 10 U.S. gpm for 2 
months. The radiation field on contact with this pipe 
when removed from the loop was ~ 20 mr/hr. 

The defect test in this loop lasted for approximately 
14 months. At the end of the test a piece of pipe which 
had been onstream the whole time was removed from 
the loop and the corrosion product film removed by 
brushing and acid etching. It was found that the relative 
amounts of fission product and corrosion product activities 
were entirely different to those listed in Table 3, viz. Co® 
accounted for 50%, and Mo% 30% of the total 
y-activity after four days. 


Decontamination after fuel element ruptures 


Fuel element ruptures occur occasionally in loop tests 
causing very high activity levels in the loop water. When 
this happens with U O, fuel elements, the activities 
released to the coolant are mainly rare gases and water 
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soluble fission products such as the kryptons, xenons, 
cesiums, iodines, etc., unless the rupture is severe enough 
to allow particulate UO, to be spread throughout the 
loop. Draining the loop followed by refilling and circu- 
lating through fresh ion exchange resins greatly reduces 
the radiation level. For example, when a UO, fuel ele- 
ment ruptured in a loop at NRX at the end of F ebruary, 
1957), the activity in the loop water rose by a factor 
of 500 and radiation fields in contact with loop piping 
were in the order of 75 roentgens (r)/hr. (normal field 
< 0.1 r/hr.). Using the above treatment for decontamina- 
tion the redietion field on the loop piping was brought 
back to normal in about a week. Simply draining the loop 
shortly after the rupture occurred reduced the field to 
about 1 r/hr. 

In another similar case of a UO, fuel element rupture, 
two other decontaminants were tried as well as water— 
namely, periodic acid and citric acid. Results indicated 
that neither of the acids was any better as a decontamin- 
ating agent than water alone. Versene has been tried at 
different pH’s as a decontaminating agent with similar 
results. 

In reactors, especially those employing heavy water, 
such draining and flushing techniques are not practical. 
Decontamination must be carried out by degassing and 
ion exchange methods. If UO, fuel is used and the oxide 
is retained in the fuel element there should be no decon- 
tamination problem. However, if UO, should be released 
to the water it will tend to deposit throughout the system 
and will create an additional clean-up problem, This is 
at present commanding considerable attention. 

In NRX itself fission product contamination of the 
heavy water occurred in July, 1955, from the rupture of 
a plutonium fuel rod. In this case the major decontamina- 
tion of the D,O was accomplished by means of ion 


exchange resins. 


The use of ion exchange resins in the maintenance of 
DO purity 

It has been found in the operation of the NRX reactor 
and of high temperature loops in the reactor that high 
water purity can be maintained most readily by passing a 
fraction (about 1%) of the main recirculating volume 
through a by-pass containing a filter and an ion exchange 
resin column. U sually a mixed bed resin is used consisting 
of a strongly acid cation exchange resin and a strongly 
basic anion exchange resin. Besides remov ing ionic im- 
purities, the resin column is itself an efficient filter. When 
the mixed bed resin is in the H-OH form, the water is 
kept at pH 6.5-7.0 and conductivity < 1 umho cn. 

In heavy water systems, the resin must be deuterized 
before being used—i.e. the water and exchangeable H and 
OH groups must be replaced by D,O, D and OD, to 
avoid degrading the D,O. This is usually done by packing 
a glass column with the resin in light water. Heavy water 
(ca. 99.5% D,O) is introduced at the bottom of the 
column and allowed to displace the H,O from the column. 
With a flow rate of 20ml/min. for a 10 cm. diameter 
column there is very little mixing at the interface of the 
D,O and H,O. When all the light water has been dis- 
placed from the column, high grade heavy water is flowed 
through (ca.99.8%D,O). When no further degradation 
of the DO coming from the top of the column. is 
observed, the flow is stopped and the resin removed and 
stored in a sealed container. The amount of 99.8% heavy 
water degraded is usually about 2 kg. per 7 kg. of resin 
deuterized. A typical resin deuterization curve is shown 
in Figure 2. 

The resin also serves to maintain the pH of the water 
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in the system at the desired level. If neutral water js 
desired, the resin is used in the H-OH form. If alkaline 
water is desired, at pH 9.5-10.5, the resin is used in 
the K-OH or Li-OH form. In the loop experiments at 
Chalk River, high pH is used almost entirely to avoid the 
buildup of crud on the fuel elements. 

The size of the resin column that is used in any 
particular system will depend on the purity required, the 
rate of production of impurities, the clean-up time re- 
quired, etc. Standard ion exchange practices are followed— 
i.e. for optimum performance the flow rate should be 2-10 
U.S. gpm per ft* of resin (minimum bed depth 24 in.), 

lon exchange resins, therefore, can serve several pur- 
poses in a pressurized heavy water power reactor: 

(1) They maintain the pH at the desired level. 

(2) They help keep fission product activity at a low 
equilibrium level. 

(3) They help keep the crud concentration in the 
water low, by adsorption of charged corrosion products, 
ionic or colloidal, and filtration of particulate crud. 

(+) They serve as decontamination agents when fuel 
ruptures occur. 


Upgrading of heavy water 

In any heavy water reactor the isotopic purity of the 
D.,O must be kept above a certain level or the reactor 
will not operate without the addition of enriched fuel. 
It is therefore necessary to have some means of upgrading 
the DO, either intermittently or continuously. 

In the NRU reactor at Chalk River it is planned to 
upgrade continuously, if this is necessary, using either a 
bubble cap plate column or a sieved plate column. The 
column will be fed from a sidestream with 99.76% D,O 
which will be returned to the reactor as 99.80% D,O. : 

The heavy water in the NRX reactor does not require 
continuous upgr ading. If accidental downgrading takes 
place the D,O is reconcentrated by electrolysis. A current 
of 1000 amperes at 3 volts is used in each of 4 cells. The 
average grade of D,O reconcentrated is 97% providing a 
product (60% of the feed) which is about 99.9% D,O. 
The gases leaving the top of the cells are recombined for 
future feed. The minimum grade of heavy water handled 
is 80% D,O. To provide passage for the large current the 
feed is made 1.5% KOH and 1.0% K.CO.. 


Conclusion 

Results from experimental work in loops have indi- 
cated that a pressurized water power reactor should be 
oper rated if possible with the coolant at high pH to 
minimize the problems created by corrosion products. 
Continuous side-stream purification with mixed bed ion 
exchange resin columns is desirable to keep the crud 
concentration in the water low, to minimize the equili- 
brium fission product activity level, to maintain the 
desired pH and to clean up water soluble activity result- 
ing from fuel ruptures. 

Some means of upgrading D,O is desirable for heavy 
water reactors for intermittent or continuous use. 

It should be emphasized that the results reported here 
have been obtained on small light water loop systems and 
these results may not be the same as will be obtained in 
a power reactor. 
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End Effect Corrections 
in Heat and Mass Transfer Studies’ 


A, I. JOHNSON’, A. HAMIELEC’, D. WARD® and 
A, GOLDING* 


Recently, considerable effort has been made to 
investigate mechanisms of heat and mass transfer to 
drops or bubbles passing through a_ continuous 
liquid medium. Many of these studies have been 
hampered by the unavailability of accurate methods 
for determining end effects. 


This paper outlines methods of plotting efficiency 
measurements for the overall transfer process so that 
both the magnitude of the end effects and the 
mechanism of the steady fall or rise period may be 
investigated. The examples are drawn from mass 
transfer with resistance inside the dispersed phase 
and from heat transfer with resistance in both phases. 
For resistance in one phase the methods are straight 
line extrapolations. 

Various mechanisms of transfer inside drops have 
been studied to test the methods of determining end 
effects. 


URING recent years several studies have been made on 
D the heat and mass transfer rates to and from drops 
or bubbles passing through continuous fluids. These 
studies have been directed at examining mechanisms for 
the transfer processes inside the drops or bubbles, or from 
these, during the steady fall or rise of the dispersed phase. 
Here “steady” infers that the transfer mechanism inside 
and outside the drop or bubble is constant w ith respect 
to the drop or bubble. 

These studies have usually been carried out by forming 
the dispersed phase, allowing it to rise or fall and then 
to coalesce (or burst in the case of bubbles) all in contact 
with the continuous phase. Different transfer processes 
are occurring during all three states of the drop or 
bubble. In order to measure the transfer during the steady 
rise or fall period, end effect corrections for the other 
periods must be made. 

There have been several mcthods of evaluating these 
end effects. This paper will review these and propose 
modified techniques for combining end effect analysis and 
examination of the transfer mechanism. As far as possible, 
an attempt has been made to obtain straight line extra- 
polation techniques. 


Previous work 


Among the first investigators reporting end effects in 
liquid-liquid extraction were Johnson and Bliss “'), Sher- 
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wood, Evans and Longcor °°), Treybal “), and West, 
Robinson, Morgenthaler, Beck, and McGregor “*), 


Sherwood, Evans, and Longcor determined a drop 
formation end effect by plotting the logarithm of the 
fraction unextracted based on the concentration of the 
solvent phase, against column height. The resulting 
straight line was extrapolated to zero column height to 
give an end effect of approximately 40% of the total 
material transferred. Using this procedure, West, Robin- 
son, Morgenthaler, Beck, and McGregor repeated the 
work of Sherwood, Evans, and Longcor and found with 
a similar plot a formation end effect of about 20%. The 
reason for this inconsistency was probably due to direction 
of transfer. : 

Licht and Conway ‘5) have indicated that a study of 
the mechanism of solute transfer should recognize three 
separate stages of extraction in the life of “each drop 
during spray tower operation. The period during which 
the drop is forming on the tip constitutes stage 1. Stage 
2 covers the period of drop fall (or rise) through the 
continuous phase. Stage 3 includes the coalescence of the 
drop and any interface transfer at the opposite end of 
the tower. They have attempted to determine formation 
and coalescence end effects separately, using a specific 
apparatus for eliminating the coalescence end effect. Licht 
and Conway plotted percent extracted for the disperse 
phase versus column height to determine end effects. 
However, for several cases encountered the plots were 
curved, but were extrapolated as straight lines. The end 
effects thus determined were probably not very accurate. 

Licht and Pansing ‘") concluded that the amount of 
extraction occurring during drop formation is sufficiently 
small so that experimental results do not detect the varia- 
tion with drop formation time. They also state that unless 
the final end effect is eliminated any extrapolation pro- 
cedure is incorrect. 

Coulson and Skinner ‘7) developed an ingenious appa- 
ratus for determining the amount of transfer occurring 
during drop formation. Drops were formed on the end 
of a capillary and then withdrawn before breaking away. 
Their technique did not measure the transfer caused by 
drop deformation as a drop is detached from the orifice. 

Pang Sheng Li ‘*) reported end effects for absorption 
from bubbles of various gases. 

Geankopolis and Hixon ‘9) measured the mass transfer 
rates in various sections of a liquid-liquid extraction 
column by sampling the continuous phase by means of a 
moving probe. They reported a large end ‘effect at the 
continuous phase entry-end. A similar study on a second 
system was reported by Geankopolis, Wells and Hawk (19), 
Gier and Hougen (11) made a similar study of spray 
columns in liquid-liquid extraction and reported large 
circulation effects caused by rising liquid drops. 
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Apparatus 

The _— used in this study is shown in Figures 
1 and 2 

Figure 1 illustrates the apparatus in which the transfer 
of n- -butyl alcohol into falling drops of water has been 
studied. The same apparatus without the linseed oil sepa- 
rating phase was used to obtain data for the transfer of 
cy clohexanol from a continuous phase of cyclohexanol 
saturated with water to falling water drops. Tt consisted 
of a glass cylinder 3 in. inside diameter and 3 ft. long 
tapered at the bottom to a stopcock through which 
coalesced disperse phase could be withdrawn. 

It was found that drops of water greater than 0.00900 
ft. in diameter would not form well in the continuous 











Four outlet 


Thermocouple glass tubes 


Nozzle 


Distilled water flow from 
. constant head tank 


Drop catcher 


Glass _ sliver 
_l> Mercury seals 


Rubber stoppers 


—Four inlet 
Thermocouple glass tubes 


Drain 







n-butyl alcohol Figure 2—Apparatus for measuring heat transfer rates 
s from or to drops, with resistance in both phases. 
saturated with 


water phase of n-butyl alcohol saturated with water and had 

to be formed in air above the continuous phase. However, 
for the water-cyclohexanol system all water drops were 
formed in the continuous phase. The single dispersed 
phase drops were formed from sharp-edged stainless steel 
nozzles 5/16”, 1/4”, 3/16” and 1/8” inside diameter, as well 
as from a number of hypodermic needles (7), Drop diam- 
eters were determined by counting the number of drops 
required to make a volume of approximately 7 cc., assum- 
ing the drops to be spherical. 


In an attempt to eliminate the coalescence end effect, 
linseed oil was carefully run in under the continuous 
> phase (n-butyl alcohol saturated with water), to act as a 
separating layer. 
| =— Addition of The amount of solute transferred to the drops was 
continuous phase measured by the refractive index change of the collected 
disperse phase. 

Figure 2 illustrates the apparatus in which the heating 
and cooling of dispersed phase falling through a con- 
tinuous phase was observed. The apparatus was essentially 
the same as that used by Korchinski “'?) except that the 
nozzle was jacketed throughout its length. It was found 
advantageous to have the drops fall on a glass point in 


Linseed oil separating 


phase the drop catcher, to decrease the time required for drops 
to coalesce and thus decrease the end effect. This column 
Coalesced dispersed had an inside diameter of 37 mm. and was 53 cm. long. 


The length of the column for drop fall was varied by 
adjusting the position of the nozzle in the column. Copper- 
constantan thermocouples were used to measure temper- 
ature. 


phase 


= . Calculation techniques 
Figure 1—Apparatus for measuring mass transfer rates 


inside drops, showing second immiscible phase for elimin- The calculation techniques discussed in this section 
ating the settling end effect. will be illustrated later by both mass and heat transfer 
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ex} riments. For brevity, the definitions and equations 
given in this section refer to mass transfer operations only. 


a) End effect definitions: The following definitions may 


ye stated: 
Ci = €, 
So te. re a 1 
‘G4 =-Cc (1) 
=f, 
a i ee oie re tie ont: 2 
aie (2) 
E vs Ce _—— C3 (3) 
bm Gates 
«=< 
Er = ——............. eee 
. y= ( 


where the various concentrations have been defined on 
Figure 3, with the exception of C* which is the equili- 
brium concentration of disperse phase reached after an 
infinitely long contact with the continuous phase. For 
any run C* is a constant. 


Period of formation, 
initial oscillation and 
drop acceleration or 

deceleration 


Steady fall 
period of drop 


Period of coalescence 
and contact of 
settled dispersed phase 
with continuous phase 





Figure 3—Definitions of transfer regions during drop fall. 


Here E,, is the fractional approach to equilibrium 
during the steady rise or fall period. It is for this period 
that one wishes to be able to test various mechanisms. 

On the other hand Ey is the fractional approach to 
equilibrium calculated by taking measurements on the 
entering and leaving dispersed phase. 

The initial end effect E;, may include several transfer 
mechanisms—transfer inside or at the nozzle and excessive 
transfer rates during the unsteady period usually occurring 
after drop release. The effect Ey, is that observed during 
coalescence and may include transfer to the settled dis- 
persed phase. Normally for experiments at constant drop 
size and formation rate E,, and E,, are independent of 
column height. 

From the definitions above, the following may be 
obtained algebraically: 

E - Er e En —_ Exe + En x Exe - Ey - Er (5) 

"1 = En — En + En X En 1-Er ' 
Where E, is the combined end effect Ey, + Ey, — Ey, 
X Eve 


In order to study the mechanism of the steady rise or 
fall of a drop it is necessary to either eliminate E,, and 
E,,, or to determine and allow for the quantity E,. 

‘If it were desirable to obtain either the fonmati or 
coalescence end effect it can be seen from Eq. (5) that 
it is necessary to eliminate one of these effects. If one 
effect is eliminated then Eq. (5) is simplified to 

3 Er — Ey 
Em = go fi eee (6) 


or 


cen. svctother . (7) 


b) Transfer mechanisms: 

During the steady fall or rise period of a drop heat or 
mass is transferred inside the drop by a transient process. 
The more important theoretical equations describing this 
transfer process are: 

1) The Newman equation “'5): negligible continuous 
phase resistance and transfer in the drop by a molecular 
diffusion process: 


_ Dn t 


Qe is 
En = lk --i- B® 


(8 
G=-¢ mr ) 


n=! 
ii) The Vermeulen equation ‘!4) is an empirical equa- 
tion closely approximating the Newman equation: 


or 


— in(1 — En?) = 


r —_ 2t. ‘ . . . 
By expanding e in the right side, this equation may 


be simplified to 
En = —Di!..... agate 


It can be shown that this fits the Newman equation well 
for E,,, values less than 0.5. 

iii) The Kronig and Brink equation ‘'5): negligible 
continuous phase resistance and circulation within the 
drop: 

co 
Em =1-15 Ble ee 
n=! 

iv) A modified Vermeulen equation has been shown 
by Korchinski “'?) to represent the Kronig and Brink 
equation well: 


where R had a value 2.25. Again, for E,, values less than 
0.5 this may be modified to the following equation similar 
to (9a) 

: RirbD! t} 

Diy, ea eateries oom odayeaees (11a) 


Ta 


v) When there is appreciable resistance in the con- 
tinuous phase the following equations by Jakob ‘'®) des- 
cribe the process: 


& 6 (an Z,:— Z, coeZy) ' — DZ,?t (12) 
ee = 2.33 3 (Z, — sin Zy cos Zn) : ra? i 


n=1 
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where Z,, are the roots of the auxiliary equation (for heat 
transfer ) 


her 
— = 1 — nrcotnr.... 
KK 
vi) For a circulating drop with appreciable outside 
resistance the following equation is given and is based on 
Korchinski’s premise that circulation streamlines act to 
effectively increase the molecular diffusivity by a factor R. 


Bed (sin Z,» — Zy cos Zn)? R DZ,? t 
Em = “27 — - l1—e oD 
a a — 6in ha cos Za) ra? 


For certain systems this factor R was found to be 2.25 by 
Korchinski 12), 


c) Plotting methods 


i) Resistance to transfer inside the drop: If there is 
resistance to transfer inside the — then Eq. (9) and 


(11) would indicate that the —In (-—E,,?) should sce 
bd : Dr? t + - 

against the group The procedure here is to assume 

£ g ra 

values of E, and convert each Ey value to E,, which is 


2 


then plotted as 1—E,,? against on semi logarithmic 


Td 
am . Dr? t : 

E,,”) against 12 On ordinary rec- 
tangular coordinates. The correct E, value is that one 
which causes that data to lie on a straight line passing 
through zero. From the slope of the rectified line the 
modifier R [such as the 2.25 value of Eq. (11)] of the 
diffusivity can be obtained directly. 


paper, or as —In(1— 


If sufficient data can be obtained at E,, values lower 
than 0.5 (for the water-cyclohexanol system in_ this 
study), Eq. (5) and (lla) can be combined to yield 

Ria D} 
= (1 — Ef) t2 es a eae C13) 
ra 
This equation indicates that E, should be plotted against 
t' on rectangular coordinate paper. The intercept yields 
E, and the slope may be used to obtain the correction 
factor R. 

ii) Resistance to transfer in both phases: Examination 
of equation 14 indicates that there is no method similar 
to the above for plotting E, data to obtain straight lines. 
If it is desired to estimate the end effect E, and the 
exponent R from a set of data, a series of values of both 
E, and R must be chosen until the data agree with the 
theoretical curves. This curve fitting technique will be 
illustrated by a sample c calculation. 

iii) Approximate method for obtaining E, for either 
case: an empirical method for wane te data on 
rectangular coordinate plots uses 


Ey’ = (1 JEr + A.... see OD 


where A is a constant with a recommended value of 0.80. 
Curved plots of Ey against t are considerably straightened 
and the intercept of the E;' against t can be used to 
estimate the total end effect: 


intercept = (1 — A 
Although limited in accuracy, this method may be useful 
and will be illustrated by a sample calculation. 
iv) Checking a suspected mechanism: Eq. (5) may be 
rearranged to give 
(1 — Er) = (1 — En) (1 — Ef)... nik eee 
or in(1 — Er) = 1n(1 — En) + In(1 — Epf).....(19) 
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Since E, is independent of time t, differentiating yields 


d in(1i — Er) 


d In (1 _ Em) 
dt dt 


This a indicates that at any time the slope of =~ 
In(I—E,) curve is the same as the slope of a In(1—F,,) 
curve. Hence, on semilog plots, comparisons of experi- 
mental data with theoretical equations may be made 
regardless of the magnitude of the combined end effects 
ES. 


’ 


Illustrations of the above methods 


Tables 1 to 4 show typical mass transfer data obtained 
by Hamielec “'7)- By using aluminum particles suspended 
in the drops the runs of A, A°, and C were shown to be 
circulating while the drops in series B were oscillating. 


TABLE ! 


Mass TRANSFER DaTA FOR THE WATER — n-BUTANOL SySTEM 
AT 25° C. SERIES A° 


Em for Er = 0.300 


Contact time Er 
(seconds) 


442 

435 | 0.199 
503 

.507 0.293 
558 

551 0.365 
590 

590 0.415 


WN 


nO 


The drop diameter was calculated to be 0.00982 ft. The correction 
factor for circulation, R = 7.69 


TABLE 2 


Mass TRANSFER DATA FOR THE WATER — n-BUTANOL SYSTEM 
AT 25°C. SERIES B 


Contact time Er 
(sec onds) 


0 0.365 
ae 0.421 
4 0.494 
J 0.553 
.6 0.598 
“0 0.650 
0 0.685 


NIOARWN ee | 


The drop diameter was calculated to be 0.0124 ft. The correction 
factor for circulation, R = 29.6 


TABLE 3 


Mass TRANSFER DATA FOR THE WATER — CYCLOHEXANOL 
SysTEM AT 25°C. SERIES A 
Contact time Er Em for Eg = 0.150 
(seconds) 
2 0.207 
3 0.195 0.060 
10. 0.270 
10. .258 0.134 
20. | 0.321 
20. .321 0.201 
30. 357 
30. . 362 | 0.247 


The drop diameter was calculated to be 0.0184 ft. The correction 
factor for circulation, R = 2.26. 
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TABLE 4 


Mass TRANSFER DATA FOR THE WATER 
System AT 25°C. SERIES C 


CyCLOHEXANOL 


Contact time Ex Em for Er = 0.110 
seconds ) 
3.0 0.188 
3.2 0.181 0.084 
9.8 | 0.222 
10.0 | 0.222 0.126 
17.8 | 0.285 
17.6 0.289 0.199 
26.4 0.320 
26.4 0.320 0.236 


The drop diameter was calculated to be 0.0224 ft. The correction 


factor for circulation, R = 2.88. 


In Figure 4 the data of Table 2 for oscillating drops 
are shown. The Ey, values were obtained by Eq. (5) from 
the tabulated E, values using the assumed values of E, 


shown. It is seen that as larger E, values are used the dun 
is straightened and a line through this data is made to 
pass through zero. The correct Ey, is estimated to be 
0.135. From the slope of the line corresponding to this EF, 
value the effective diffusivity was found to be 28.3 x 10“ 
or 29.6 times the molecular diffusivity. 


ft?/sec. 


10 





0.8 


0.6 


— In (I—E?) 


Assumed E, values 


0.2 
a 
S €E,= 0.100 
& E-= 0.130 
True E-= 0.135 
% | 2 3 4 5 6 € 


Contact time sec. 


Figure 4—Observed and corrected mass transfer data for 
the water-n-butanol systems, Series B (Table 2) employing 
an end effect technique indicated by equation 9. 


In Figure 5 the data of Table 3 are plotted as E, 
against t’* in accordance with E q- (11a). The straight tho 
has been extrapolated back to zero time to yield an E, 


The Canadian Journal of Chemical Engineering, October, 1958 





0 | 2 3 4 5 6 
+2 - seconds? 


Figure 5—Observed mass transfer data for the water- 
cyclohexanol system, Series A (Table 3) exhibiting an 
extrapolation technique indicated by Equation 15. 


value of 0.150. From the slope of this line the R value was 
found to be 2.26. 


The data of Tables 1 and 4 may be plotted on graphs 
similar to Figures 4 and 5 to obtain ‘the R values indicated. 


1.0 





e 6©Ey data 
Ew for Ee s 
Ew for Er s 


0.30 


0.45 
08 


06 


04 


Theoretical curve for 
circulating drop 


0.2 
Theoretical curve for 


stagnant drop 


Contact time sec. 


Figure 6—Observed and corrected heat transfer data for 
the bromobenzene-aqueous glycerine solution system (Table 
5) exhibiting a curve fitting technique. 
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TABLE 5 
HEAT TRANSFER DATA FOR BROMOBENZENE DROPS BEING 
CooLeD BY FALLING THROUGH GLYCEROL SOLUTIONS 







Ey Contact time 


(seconds) 
me 


0.594 .90 
40 


6: | 3.00 
0.691 3.80 
0.698 | 4.00 
0.722 | 4.50 
0.757 | 4:68 


The correction 


The drop diameter was calculated to be 0.635cm. 
factor for circulation, R = 2.25 


In Figure 6 the theoretical equation of Jakob for the 
cooling of non circulating drops [Eq. (12)] and the modi- 
fied equation [Eq. (14)] using R = 2.25 have been plotted. 
The E, data of Table 5 and the corresponding Ey data 
obtained using the assumed E, values indicated are plotted 
as points. It is seen that as increasing E, values are used 
the data points are moved downward. At E, = 0.30 the 
data became superimposed on the Eq. (14) line. Within 
the accuracy of the data theoretical lines for other R 
values of Eq. (14) did not fit the data regardless of E, 
values used. 


1.0 


08 


06 


04 


0.2 


e@ €, data 


© Ey, data at As0.90 


Contact time sec, 
Figure 7—Observed heat transfer data (Table 5) indi- 


cating an extrapolation technique for estimating the order 
of magnitude of end effects present. 


Finally, Figure 8 illustrates the comparison of data 
with theoretical equations. In this diagram it is easily seen 
that the experimental curve is parallel to the predicted 
curve of the Kronig-Brink equation. This indicates that 
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Oo Experimental data 


— Kronig-Brink equation 















8 12 16 20 


Time - seconds 


Figure 8—Observed mass transfer data for the water- 
cyclohexanol system, Series A (Table 3), exhibiting a 
mechanism test. 


24s 28 32.36 


the Kronig-Brink equation holds for the steady fall period 
over the range of the experimental data. 


Contact time for mass and heat transfer studies 


It should be noted that the contact time for the steady 
fall or rise of a drop is less than the time required from 
drop break off to coalescence. The time to be used in the 
steady transfer equations should be 


t= t! — At 


where At is the time of the unsteady fall or rise period. 
Eq. (9) should be modified to 


Dr?t! Dr? At 
— in(i — E,?) = — — —— 
ra? ra? 
agli, lay ; ‘ Dr? 
This indicates that a negative intercept of value “_ 
d 


should be expected when the data are plotted in the 
manner shown in Figure 4. However in the work en- 
countered in this research to date, straight lines passing 
through the origin have represented the data well. It must 
be concluded that the time of the unsteady fall period is 
small compared to the time of the steady fall period. 


Conclusion 


This paper has attempted to define and to indicate 
simple calculation methods for obtaining the end effects 
for mass or heat transfer to drops or bubbles formed 
passing through a continuous medium from a nozzle to a 
settling area. Such end effects must be known before 
mechanisms of transfer during steady rise or fall periods 
can be investigated. Although the calculations are simple 
the present authors believe that they have been overlooked 
by most previous workers in this field. 


Nomenclature 
A = empirical constant 
B= coefficient 
C = concentration of dispersed phase 
D = molecular diffusivity 
De = effective diffusivity 
E =a dimensionless number whose value is given by (Em < 
E < Er) 
h = adimensionless number = — in (1 — Er) 
h = film transfer coefficient 
K = effective thermal conductivity inside drop. 
R = correction for circulation or oscillation of drop phase = 
Po. ; 
, dimensionless 
D 
r = drop radius 
At = time from point of drop release from the nozzle to the 
point when transfer mechanism of steady fall or rise 
period holds 
= contact time of drop for steady fall or rise period 
t! = time taken by drop to strike the coalesced dispersed phase 


from the point of drops release from the nozzle 
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Greck Letters 


Lb eigenvalue 

Subscripts 

1 at column inlet 

2 = beginning of the transfer mechanism of the steady fall or 
rise period 

3 = point at which falling or rising drop strikes the coalesced 
dispersed phase 

4 = at column outlet 

T = overall 

c = continuous phase 

d = drop or dispersed phase 

m = steady fall or rise period 

n = term number in series sum 

f; = end effect at entry end for dispersed phase 

f. = end effect at exit end for dispersed phase 

Superscripts 

* = at equilibrium with respect to transfer of solute to or from 


dispersed phase 
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LETTER TO THE EDITOR 


New York, N.Y. 





Editor, C.J.Ch.E.: 
RE: Viscous Flow in Multiparticle Systems 


This abstract reviews a series of theoretical stud- 
ies which has resulted in achieving a better under- 
standing of such phenomena as resistance to flow of 
fluids through beds, sedimentation of small particles 
and viscosity of suspensions. There are many systems 
both in nature and technology involving these types 
of motion. Examples of commercial applications 
include the separation of particles by centrifuging 
and sedimentation, filtration, flow of gases through 
vessels containing packed solids (kilns, gas producers 
and blast furnaces), and motion of particles and 
gases in catalyst contracting vessels. 


In order to construct tractable mathematical models 
for the complicated flow systems involved, it is necessary 
to resort to some simplifications. Here, the inertia terms 
are omitted from the well known Stokes-Navier equations 
of motion, resulting in the so-called creeping motion 
equations. The operations noted above are often carried 
out in the Reynolds number range (based on particle 
diameter) of five or less and consequently are amenable 
to this assumption. The other, important simplification is 
that the systems consist of uniform sized, spherical par- 
ticles. The justification for these assumptions rests in part 
on the conciseness of treatment thus made possible and in 
part on the good agreement of results with available 
experimental data. 

Essentially two techniques have been employ ed, 
namely the ‘method of reflections and the unit cell. In 
the method of reflections the equations of motion are 
satisfied successively in each of the separate bounding 
surfaces including the container walls of a suspension. If 
the system is sufficiently dilute, rapid convergence is 
possible. The unit cell technique i involves the concept that 
an assemblage can be divided into a number of identical 
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cells, one sphere usually occupying each cell. The 
boundary value problem is thus reduced to a consideration 
of the behavior of a single sphere and _ its bounding 
envelope. This technique applies best where symmetry is 
more or less complete in the particle assemblage and thus 
is of greatest applicability in concentrated assemblages 
where the effect of container walls can be neglected. 

An approximate theory for the behavior of dilute 
multiparticle systems suspended i in a viscous fluid has been 
developed using the reflection technique “). This theory 
is based on a rigorous treatment for the case of a single 
sphere free to occupy any position, both axial“) or off- 
axis (*) in a cylindrical tube. The results obtained include 
estimates of particle velocity and distribution as well as 
pressure drop due to fluid flow for several systems. They 
are in agreement with the qualitative data av -ailable. 

In the case of the unit cell technique a novel approach 
has been developed on the basis that each of the particles 
is associated with an envelope of surrounding fluid which 
in turn is bounded by a “free surface”. This model enables 
closed solutions of the hydrodynamic equations to be 
obtained. From these results rate of sedimentation or 
alternatively pressure drop can be predicted as a function 
of fractional void volume“). In the case of shearing 
motion similar solutions relate relative viscosity of a 
suspension to solids concentration“). Agreement with 
available data is very good without the use of any 
empirical constants for both resistance to flow and vis- 
cosity over a wide range of concentration. 

In short, these studies have resulted in useful models 
for extension to more complicated cases, aside from their 
possible immediate application in practical correlations 
based on experimental data. 
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